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ABSTRACT 


The omni-present low frequency wave motion (30-300 sec) contains a substantial 
fraction of the total wave energy inside the surf zone. A more complete description of 
nearshore wave processes considers incident short period wind waves superposed on, 
and interacting with, long standing waves. The wind waves are modulated in 
amplitude, wavenumber and direction due to relatively slowly varying depth changes 
caused by the long waves. The energy in the wind wave band is enhanced by side band 
growth at the sum and difference frequencies of short and long waves (order 15% at 
the shoreline). The modulation is identified in the analysis of field data as a positive 
correlation between the long waves and the wind wave envelope near the shoreline. 
Considering oblique incident waves, a _ steady longshore current showing a 
non-vanishing current at the shoreline is found as a result of the non-linear interaction 
between monochromatic incident and infragravity waves. An analytical solution 
describing the unsteadiness of the longshore current is developed. Applying the derived 
longshore current solution, longshore sediment transport 1s reformulated to include the 


infragravity waves, giving improved comparisons with field measurements. 
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I. INTRODUCTION 


Sca waves are generated by wind, and as they propagate towards the shore they 
deform due to shoaling and may refract if they are oblique to the bottom contours or if 
they mect a current. [The waves grow in height until they beGoOme Mista wle amd: orean. 
dissipating their energy. The processes of breaking waves on a sloping bottom make 
the surf Zone an extremely dynamic area where different observable phenomenac take 
place, such as longshore current, sediment transport and wave runup (setup and 
swash). Due to this complexity, the proper modeling of wave motion in the nearshore 
zone has been the goal of many investigators, and yet much remains to be done before 
a satisfactory understanding is achieved. 

Miiessuit Zoncecamebe Charactenzed dynamically by three regions: the outer,the 
iM@ereana the tunup resion (Svendson et al., 1978), (see Fig 1.1). The outer region 
follows the breaking point and ts dominated by incident breaking waves. The breaking 
waves are of either the plunging or spilling type, and are characterized by a rapid 
transition of wave shape with a horizontal surface roller. In the inner region, the wave 
forms break down into smaller seales of random and turbulent nature and eventually 
become similar to moving bores, or hydraulic jumps. The run-up region is 
characterized by wave setup and swash oscillation running up and down the beach 
face. In all regions, the wave height and the water motion are strongly locally 
controlled by the depth. The depth controlled waves are referred to as “saturated.” 

In general, surf zone dynamics are found to differ dramatically between reflective 
(steep) beaches and dissipative (gentle slope) beaches (Bradshaw, 1980). Typical 
features of the reflective beaches are high reflection of incident waves, collapsing or 
surging breakers on the beach face, and subharmonic resonances (c.g. Guza and Davis, 
ei za and Inman, 1975; Eluntley and Bowen, 19756, 1978; Wright et al., 1979). 
Dissipative beaches are dominated by spilling breakers, pronounced low frequency 
oscillations and the presence of multiple parallel longshore bars (¢.g. Huntley and 
pereeneulo7 oa Short, 1975; Sasaki and Tlorikawa, 1975; [folman et al., 1978; Symonds 
et al., 1982). This study will focus on the dynamics of dissipative beaches in which the 
inner region 1s assumed to be a transition region between regions dominated by 


incident breaking waves and low frequency swash. 
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Figure 1.1 Surf zone regions and coordinate system. 
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In the past decade, considerable attention in the literature of beach and 
nearshore dynamics has been devoted to the kinematics of low frequency waves. These 
waves, Which have periods ranging from 30-300 sec, are referred to as “infragravity 
waves. The observations show that the dominant wave motion closer to the shore is 
not normally at the incident wave frequency band. It is also found that the low . 
frequency wave energy in the very nearshore region contains a substantial fraction of 
the total wave energy (Inman, 1968 a,b; Suhayda, 1974; Goda, 1975; Huntley and 
Bowen, 1975b; Sasaki and HTlorikawa, 1975,1978; Iluntley, 1976; Wright et al., 
1978,1979,1982; Bradshaw 1980; Holman, 1981; Huntley et al., 1981; Guza and 
Thornton, 1982,19$5a,1985b; and others), (see Fig 1.2 for example). The valley in the 
energy spectra, Where the energy level reaches minimum, suggests an empirical cutolf 
frequency around 0.03 Hz that separates the two wave bands. Therefore, a complete 
description of the wave field in the surf zone should include two relatively separated 
bands of wave frequencies : (1) incident sea-swell waves (short period waves) with 
periods of 1-30 sec, and (2) infragravity waves (long period waves) with periods of 
30-300 sec, which no longer can be neglected inside the surf zone. 

Prior to the 1950's, the incident wave characteristics were the main concern when 
measuring and analyzing waves in the nearshore zone. Munk (1949) and Tucker (1950) 
were the first to point out the existence and significance of low frequency motion and 
called the phenomenon “surf beat.” More recently, with renewed interest 1n nearshore 
processes and with major improvements in the instrumentation, 1t has become apparent 
that the low frequency wave energy can be as important as the incident wave energy 
when approaching the shoreline. Another motivation for the study of the infragravity 
waves is the coupling between the infragravity length scales (order of hundreds of 
meters) and the length scales associated with most of the rhythmic beach features such 
as sand bars, beach cusps and crescentic bars. Many imvestigators have related the 
beach morphology changes to the existence of low frequency motions by which a 
mechanism of forming and maintaining such topographic features is provided (Bowen 
and Inman, 1971; Guza and Inman, 1975; Bowen, 1980; Holman and Bowen, 1982: 
Bowen and [luntley, 1984; and others). 

In analyzing surf zone data, there are often anomalous patterns of water motion 
closer to the shorelme when compared with present short wave models. [I’or example, 
the spatial distribution of wave height measured in both the laboratory and ficld often 


shows unexpected wave height at the shoreline (e.g. Suhayda and Pettigrew, 1977: 
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Battjes and Stive, 1985), (sce for exaimple Fig 1.3). Bowen, et al. (1968), in measuring 
the setup, found a residual wave height always present at the shoreline and modified 
their sctup formula at the beach face to fit the data. Suhayda and Pettigrew (1977) 
measured the average wave velocity across the nearshore zone and found it to decrease 
shoreward until the sea-swash limit of the run-up zone was reached, and then to 
increase onshore, supporting Waddell’s measurement (1973). Thus, the wave height 
and celerity do not vanish at the vicinity of the shoreline, invalidating the widely 
apphed short wave saturation assumption. 

The assumptions of wave saturation and steady state conditions for both 
monochromatic and random waves, results in a zero current at the shoreline. In 
addition, the steady longshore current assumption is often unrealistic, especially with 
the presence of irregular oscillating waves (Wood and Meadows, 1975). This is often 
seen in field data that show longshore velocity oscillations at longer wave periods than 
for incident waves (c.g. Holman and Bowen, 1984; Guza and Thornton, 1985a; 
Oltman-Shay and Guza, 1986). Using an empirical orthogonal function (EOF) analysis 
on 64 data sets covering a wide variety of wave conditions, Guza et al. (1986) obtained 
objectively best fit to the longshore current distribution. The analysis showed the 
“classical parabolic” shape for the longshore current distribution (Bowen, 1969) over 
most of the surf zone, but with significant residual velocity at the shoreline (Fig. 1.4). 
In the analysis of the same data, Wu et al., (1985) found a poor agreement between 
longshore current theory with velocity going to zero at the shoreline and data closer 
to the shoreline. All present models predict Zero current at the shoreline, and no model 
is available to explain the unsteadiness of longshore current. 

Improving the longshore current description 1s important because it provides the 
driving means of transporting sediment alongshore (Longuct-lliggins, 1972). The 
available longshore sediment transport models (e.g. Byker, 1971; Thornton, 1973; 
Komar, 1977; Sawaragi and Deguchi, 1979; Bailard, 1981), all show a longshore 
transport maximum closer to the breakerline, which agrees with the longshore current 
distribution. The shortcoming of these models is their failure to predict anv moving 
sand at the shoreline regardless of the wave conditions. Zenkovich (1960) measured 
the suspended sediment transport distribution across the surf zone using fluorescent 
sand tracers. Sand transport was found greater over bars where energy dissipation is a 
maximum due to wave breaking and also at the shoreline, (see Fig 1.5). Mlcasuring 


longshore sand transport distribution using sand tracers, Kraus et al. (1981) found a 
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distinct bimodal distribution having niaxima at the breakerline and nezr the shoreline. 
The same results were obtained by Downing (1984) when measuring the suspended 
longshore sediment distribution. Abdelrahman (1983) compared longshore sediment 
transport at Santa Barbara beach with his predicted longshore sediment transport 
model based on the widely accepted energetic concept by Bagnold (1963). The model 
did not predict the large transport in the vicinity of the shoreline since theoretical short 
wave description gives zero amplitude at the shoreline. It was suggested a primary 
cause of these differences was the presence of infragravity waves, which have not been 
taken into consideration before. 

In recognition of the importance of the infragravity waves in the surf zone, an 
obvious next improvement in the description of waves kinematics is to combine the 
infragravity Waves with the incident sea-swell. As a result, the incident waves would be 
modulated in amplitude, wavenumber and direction. The purpose of this dissertation is 
to improve the understanding of nearshore processes by generating and applying such a 
composite wave description. 

As will be seen, the interaction of the short waves with the long waves results in 
energy not only at their respective frequencies, but energy transferred to side bands of 
the incident waves at the sum and difference frequencies. It is the energy and 
momentum at the side bands that is primarily responsible for the inner surf zone 
dynamics. The dynamical changes are demonstrated by the successful treatment of 
using the first order wave description in the momentum flux equation to obtain a 
second order dynamical effect. Ilowever, inclusion of infragravity waves does not 
mean introducing more independent parameters since the infragravity wave energy 1s 
assumed related to the incident wave field. 

In the following chapter, a description for the short and long period waves 1s 
presented. The short period waves are described by linear (Airy) theory, while the long 
period waves on a sloping bottom are described by either cross-shore standing waves 
(2-D) or edge waves (3-D). In chapter 3, the changes in the amplitude and 
wavenumber of the short waves superposed on long waves are discussed. An analytical 
expression describing the modulated surface elevation of the short waves is developed 
along with their energy and momentum. Changes in the wave kinematics are 
considered with the long period waves treated as a slowly varying current. In chapter 
4, narrow banded wave data are analyzed to examine the hypothesis of short wave 


modulation by the long waves. The dynamical effects of considering the infragravity 


17 


0.4 





0.2 
, DISTANCE OFFSHORE t 
‘oO . Ore 0.4 
a 0 
E ort mms 
I> 
BeQee 
“apa 


—— first eignfunction 


~---- second eignfunction 


figure 1.4 Longshore current distribution using empirical orthogonal function 
analysis. The first crgnftnction contains most of the varrance 
(alter GiuZzaer dilaton 


18 


Waves inside the surf zone are investigated. Analytical models for calculating wave 
setup, longshore current and sediment transport are developed in chapter 5. An 
analytical solution to the unsteady longshore current is derived. Model validation and 
comparison with other models are included in chapter 5. The conclusions of this study 


icepresented in Chapter 0: 
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tia VS IN TIE NEARSITORE ZONE 


A wave record of incoming waves usually contains short period waves (1-30 sec) 
and long waves (30-300 sec). To study first order modulations of short waves 
superimposed on long waves, the principle of superposition is applied. ‘Therefore, the 
non-linear wave theories are excluded in this study. A chosen Cartesian coordinate 
system 1s illustrated in Pig. 1.1 where the positive X-axis is directed offshore, thie 
positive Z-axis is vertically upward and Y is alongshore. This coordinate system ts 
convenient for desenbing the offshore decay of the standing long waves. Linear 
descriptions of both sea-swell and infragravity waves are discussed separately. Short 
waves are derived from lincar wave theory over a flat bottom, then using shallow water 
approximation to describe short waves in the nearshore zone. Infragravity waves are 


derived from shallow water equations on a sloping bottom. 


A. SHORT PERIOD WAVES 

In general, Waves in a viscous fluid propagate over irregular topography of 
varying permeability. In most cases, the main body of the fluid motion is reasonably 
assuined irrotational sinee the viscous effects are usually lumted to a thin boundary 
layer near the surface and the bottom. Surface waves are considered a boundary value 
problem. 

Linear wave theory (small amplitude waves) 1s considered a first order 
approximation to the theoretical descriplion of wave behaviour. It 1s assumed that the 
wave amplitude , a, 1s very small compared to both the wavelength Land the local 
water depth h, (a<< LL and a<< h) The fluid is assumed incompressible, 
homogenous, inviscid and also irrotational, by which a velocity potential @ should exist 


to represent the field of the flow: 


cy 
i] 
1 


V9 (2.1) 


where u is the velocity vector. [quation 2.1 satisfies the continuity equation V .u=0, 


by which the governing differential equation (Laplace equation) 1s obtained: 
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With the proper kinematical and dynamical free surface and rigid Mat bottom boundary 
conditions, the method of separation of variables is used to solve equation 2.2 Then, 


the velocity potential 1s given by 


a = =e cosh GOs oot ai) (2.3) 


where a is the wave amplitude, 0 = (27/T) is the radial frequency, and k = (27/T.) 1s 


the wavenumber. ‘The surface elevation 1) 1s periodic in both space and tine. 


nh "aera cos (k, x; + wt) , ny ele ae (2.4) 
Where 1 = 1[,2 refers to conditions in the X and Y-directions. The description of the 


Wave mouon 1s completed by the dispersion relationship, 


w = gk tanhkh (2.5) 


Which arises from the dynamical free surface boundary conditions. The horizontal 


velovity under the wave ts given by 


7 Cosh wis( hz) - 
(US Beers jem Cos =(kse + it) (2.6) 


and the vertical velocity 


Sinn k(h+z . 2 
es au) sane te) SinUese + int) (2.7) 


Swell waves propagating from decp water are nearly sinusoidal with tong low 
crests. AS Waves approach shallow water, they start to feel the bottom and deform at 
relative depth h/L < 1/2. The wavelength decreases while the wave height generally 
Increases, and the wave period remains invariant. The hyperbolic functions have 
convenient deep and shallow water asymptotes by which an approximate solution can 
be obtained. In shallow water, the lincar gravity waves arc non-dispersive since the 


dispersion relationship reduces to 


and the horizontal and vertical velocities simplify to 


fee a yay cos (kx + wt) i) 


aw (1 +e) Sin (kx + wt) 
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(2.10) 


On dissipative beaches, long crested swell waves break when the wave steepness 
increases and the velocity of the water particles at the crest exceeds the phase speed of 
the wave form. In the inner surf zonc, bores retain their relatively long crested form 
and progressively decrease in height as they advance onshore. Due to turbulence and 
non-lincaritics associated with the breaking process, the motion of the water particles 


after breaking can no longer be deseritbed analytically. As a first approximation, the 
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wave or bore height, 11, at any location inside the surf zone 1s limited by the local 


depth according to the saturation assumption, 1.¢. 


she ile (2.11) 


where y is a constant of order unity. This relation holds for spilling breakers, the most 
common breaker type on natural beaches. The breaker index y 1s important in 
describing the breaking waves used in longshore current and sediment transport 
models. This linear relationship is expected on dimensronal grounds (Longucet-l{figgins, 
1972) and gives a reasonable description of laboratory results for monochromatic 
waves (c.g. Galvin and Eagleson,!965) and field measurements during conditions of 
both broad and narrow spectral distributions of wave energy (Thornton and Guza, 
1032) ieee): 

Very close to the shoreline in the run-up region, the fluctuating component of the 
Wave runup about the mean water level 1s known as “swash oscillation.” Miche (1951) 
hypothesized that monochromatic waves within the surf zone may be composed of 
both a progressive component, Which is dissipated onshore and has zero shoreline 
amplitude, and a standing component, which has its maximum amplitude at the 
shoreline. Many laboratory studies (Moracs, 1970; Battjes, 1974; Guza and Bowen, 
1976, Van Dorn, 1976) report that the shoreline swash amplitude a, 1s limited by the 


swash parameter €., 





where tan f is the beach slope and &, varies from one experiment to the other in the 
range from | to 3. Guza and Vhormiton (1982) attriletited this variation to tie 
ill-defined nature of the backwash (run-down). Hunt (1959) related the total runup RY 


, defined as the sum of setup and half the swash amplitude, to the surf similarity co 


with an empirical constant M: 
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Frgure 2.1 Wave heights inside the surf zone defining breaker index 
and showing saturatron (after Thornton and Guza, 1982 
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RY /Hgig = M & = M tan B/ V (Hy / Loy (2.13) 


where H.., is the significant wave height and the subscript 0 denotes the deep water 


g 
conditions. 

Huntley et al. (1977) showed a swash “saturation” condition at incident wave 
frequency, 1.¢e. with increasing wave height, the steady set-up will increase but the 


swash amplitude will not. Guza and Thornton (1982) found a «9 


spectral decay at 
incident wave frequencies, indicating that energy levels are independent of incident 
wave height, although Huntley et al. (1977) have shown 4 dependency. The swash 
saturation condition is found to be true for both monochromatic and random swash 
fluctuations (Guza et al., 1984). Swash at low frequency is found to be unsaturated, 
l.e. the swash amplitude increases with increasing incident wave height. Most of the 
low frequency energy is found to have infragravity wave periods. Thus swash spectra 
show a saturated region at incident wave frequency and an unsaturated region at low 
frequencies that dominates the whole spectrum. 

To evaluate the amount of runup, swash studies have been mainly concerned 
with measuring the swash amplitude at the shoreline without putting much emphasis 
on the generation and types of the domunant tow frequency motions. It is of interest 
now to investigate the forms of this low frequency water motion that give maximum 


amplitude at the shoreline. 


B. LONG PERIOD WAVES 
The long period waves presented in this work are limited to waves having a time 


scale of 30-300 seconds where gravity 1s the primary restoring force. 


1. Infragravity Wave Models 
The principal wave types that may contribute to the low frequency motion on 
beaches can be summarized by their mathematical formulation as follows: 


(i) 2-D model : the wave motion is mainly in the cross-shore direction and 1s 
known as surf beat or “leaky modes. 


(1)  3-D model : the wave motion has longshore variation and is known as edge 
Waves or “trapped modes”. 
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The mathematical forms of both waves are derived from the linear shallow 


Water Wave equation on a sloping bottom 





2 
moe 30 3 3b ia (2.14) 
ae? gizy(h en - ay th are a5 6 


In the 2-ID model , the wave energy is reflected offshore and the solution 
obtained 1s a standing wave expressed in terms of the zero order Bessel function of the 
first kind, Jp , (Lamb, 1932; Fredrichs, 1948) 


g(x,t) = -=2 5, (x)sin wt (2.15) 


(7 . . 
where % = (40y? x/ g tan gt “and ais the amplitude at the shoreline. It is noted that 
the Bessel function of the second kind, nae is excluded since it is unbounded as waves 


approach the shoreline. The surface elevation and the associated orbital velocities are 


given by 


s 
n(x,t) = a J, (x) cos wt (2.16) 
Meet) = avg/h J, &X) Sinwt ssl) 
(SS ers Jy (4) Sasa Coal Se) 


In the 3-1) models, the edge waves are free waves that propagate along the 


Soa iieiicimeenctyy (tapped Onshore by relraction, 1.¢. they do not radiate energy 


offshore. Stokes (1846) was the first to provide an analytical solution to describe edge 
Waves using a small amplitude wave assumption. Eckart (1951) solved equation 2.14 
analytically in terms of Laguerre polynomials L,, (2k, x). The solution of progressive 


cdlge Waves Is 


= hex 
gee a9. Y 1, (2 k x)cos(kyy - wt) (2.19) 


where ky is the longshore wavenumber, and n is an integer number which represents 
the edge wave mode number. Each discrete mode contains finite energy, and the n-th 
mode edge wave will have n= zero-crossings after which the amplitude decays 
exponentially. [idge waves obtained by [Eckart (1951) must satisfy the dispersion 


relationship, 


Depeche, (ise aaiceaas p (220) 


which requires (2n-+ 1 )tanB << 1 to insure that the solution achieves its limiting value 
as ky X —p~ 88 while still in shallow water. 

Ursell (1952) was able to obtain a set of “exact” edge wave solutions on a 
sloping bottom using the small amplitude wave theory without recourse to the shallow 


Water approximation. [lis solution satisfies a slightly different dispersion relationship 
W HG oe Qa 


and the solution hes within the range ((2n+ 1B ) S 2/2 to insure a cutoff mode for a 
given beach slope. Ursell (1952) showed that Stokes (1846) edge wave solution is only 
the zero mode of his solution. For gentle beach slopes and low mode numbers, both 
dispersion relationships 2.20 and 2.21 are in agreement, suggesting that the shallow 


Waler approximation ts appropriate under these conditions, (Guza and Davis, 1974). It 


is noted that the approximate Eckart (1951) solution is more commonly used in 
comparisons with ficld data due to its simpler mathematical formulation. In the field 
data analysis, a difficulty in discriminating between Ieaky and edge wave modes for 
small offshore distances arises because of the similarity of their cross-shore structures. 
Also, both waves have maximum amplitude at the shoreline. Fig. 2.2 shows such 
sinularities between higher mode edge waves and reflected waves. 

In summary, the low frequency motions may be classified as edge waves, leaky 
modes, forced waves or a combination as shown in Tig 2.3, (after Symonds, 1982). A 


2 


complete set of free wave solutions of discrete edge waves occurs for @* < gk, while a 


continuum of Icaky mode exists for 0? > gk It is noted that the forced waves may 


vy" 
lie between the edge wave modes, but they do not satisfy the dispersion relationship 
and are weaker than the free modes since edge waves may grow resononantly as they 


are excited. 


2. Generation of Infragravity waves 

Munk (1949) and Tucker (1950) suggested that the long period waves may be 
caused by an excess of mass carried forward by groups of high swell. The swell waves 
are assumed to be destroyed on the beach, and the excess mass transport 1s reflected 
back as a free long wave. They tested this hypothesis by comparing offshore pressure 
records, delayed by the calculated travel time of the incident wave envelope. Tucker 
(1950) observed a negative maximum correlation between incident wind wave groups 
and the rnfragravity waves at a time lag equivalent to the travel time of a forced wave 
moving with the wind wave group velocity and traveling back upon reflection as a free 
wave. 

Longuct-Higgins and Stewart (1962) showed that the variations 1n_ the 
radiation stress duc to incident wave grouping drive a forced wave such that a 
depression of the mean water level (MWL) will occur under incident high waves, and a 
corresponding rise in MWL occurs under low waves. This forced oscillation of second 
order is associated with groups until the incident waves break onshore and release a 
free wave offshore, which theoretically supports the above observation. 

Gallagher (1971) suggested edge wave excitation through a non-hnear 
Interaction whereby two incident waves can transfer energy into an edge wave if the 
difference frequency and longshore wavenumber of the two interacting wind waves 
satisfy the edge wave dispersion relationship, equation 2.21. This resonance has been 


observed in the laboratory (Bowen and Guza, 1978). 
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Iigure 2.2. Modes of edge waves and reflected waves, @~ is the velocity potential. 
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Incident waves may transfer energy to either leaky or trapped modes through 
many possible interactions. One of these is given by Guza and Davis (1974) where 
energy is transferred into edge wave modes through a non-linear interaction involving a 
normally incident standing wave on a sloping bottom. Bowen and Guza (1978) showed 
energy transfered into leaky modes (o> = $k. ) by edge waves. 

Another mechanism for low frequency Wave generation is suggested by 
collision between the stronger component of backwash and the slioreward moving 
bores that results in roll waves which last for few seconds. On analyzing swash data, 
Waddell (1973) found that the collision between uprush and backrush plays an 
important role in reducing the extent of swash runup waves. Furthermore, he 
suggested that water percolation into the beach serves as a low-pass filter affecting the 
swash frequency. Bradshaw (1980) suggested bore-bore capture as a mechanisin for 
long wave generation in the swash zone. Mase and Iwagaki (1984) showed a 
considerable shoreward decrease in the ratio between the number of runup Waves to 
the incident waves. 

Symonds et al. (1982) developed a 2-D model for generating long waves which 
is forced by a time-varying break-point caused by the groupness of the incident wave 
field. The generated free long wave at the group frequency has a standing wave 
structure shoreward of the breakpoint and an outgoing progressive wave structure 
seaward of the breakpoint. 

In the above survey, both short and long wave models are described separately 
showing different characteristics. It 1s of interest now to combine the two waves and 
study the changes in the incident wave (short wave) field due to the presence of longer 


period waves. 
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Ht. SHORT WAVE MODULATION DUE TO LONG WAVES 


When short period waves ride on much longer waves, they are modulated in 
amplitude, wavelength and direction. Their amplitude tends to increase at the crests of 
the long waves and decrease at the troughs. ‘Vhe changes in wavenumber and direction 
depend mainly on kinematical considerations. In this chapter, the changes in 
wavelength, direction and amplitude are investigated. A development for the 


corresponding changes in energy and momentum flux then follows. 


A. CHANGES JN THE WAVE AMPLITUDE 
The amphtude modulation of a short progressive wave riding on a much longer 
wave wus pointed out by Unna (1942, 1947). Unna (1947) showed a contraction of the 
wavelength and an increase in the amplitude of short waves at the long wave crests. 
Longuet-Iliggins and Stewart (1960) used a perturbation analysis to examine the 
non-lincar interaction between short and long progressive waves over a horizontal 


bottom. They described the surface elevation 


=) 
iI 


a. sin. + a, sin Vo (3.1) 


where w = (kx- @ t+ @ ), 1s the phase function and 8 the phase shift. The subscripts 
s and @ refer to the short and long wave respectively. Using Stokes’ method of 
approximation to the second order, they derived a general form for the modulated 


surface elevation 


ay @= a (1 #P)sin p. + a, QCOS es (3.2) 


where P = a, Kp sin My and Q = -ayk, cos W . Using the first term in the right hand 
side of equauon 3.2, Longuet-Iiggins and Stewart expressed the modulated short wave 


amplitude as 


a0 


i] 3 1 , th | 
a’ = a {1 + a,k, (7 cothk h + | tanhk,h)sin yo} (a3) 
and a modulation in the short wave wavenumber, k’ , as 


ae ae lal cothk ,h siny,} (3.4) 


u 
In shallow water, equations 3.3 and 3.4 reduce to 


3 ; 
a' = a_il + 7 py Sin ,! (355) 


I 
x 
= 
i— 
Bs 


k! — siny,} (3.6) 


These amplitude and wavenumber modulations are explained by the work done by the 
long waves aguinst the momentum of the short waves. ‘This work is converted into 
short wave cnergy and therefore produces a steepening of the short wave at the crests 
of the long waves, 1e., the energy is redistributed along the wavelength of the long 
wave. Longuet-Iliggins and Stewart (1960) came up with equations similar to 3.3 and 
3.4 for a progressive short wave superposed on a long period standing wave in which a 


was described as uy ; 


I. Determination of the Modulated Wave Profile 
It 1s of interest now to show that the refractive changes in amplitude and 
wavenumber due to the non-linear Interacuion may result in a modulated short wave 
expressed as an infimite sum of components at different frequencies. Then, the 
properties of the modulated wave are investigated to provide a better understanding of 
the modulation process. The procedure is to Iet the modulated short wave, designated 


by primes, propagate in the same wave direction and satisfy 


far a’ sim (kee Gee (3.7) 
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Substituting equations 3.3 and 3.4 into equation 3.7, and letting @ = Oo. , ylelds 


ion a’ sin (y. + m, sin, ) (3.8) 


Where m, = ag kp k, xX coth Kp h 1s defined as a non-dimensional wavenumber 
(wavenumber modulation factor). Equation 3.8 is a general form of amplitude and 
wavenumber modulation. [equation 3.8 reduces to the simple form H = a sin Wy 


When my, << 1. Expanding equation 3.8 using trigonometric identities, gives 


am) 6=~Ca {sin vy, cos(m, siny,) + cos p_ sin (m, sin.) } (3.9) 


Using Vourter series expansions to express the sinusoidal arguments (Abramowitz and 


Stegum, 1964), gives 


n' = a' ) J (m)sin(y. + ny,) (3.10) 


In the above analysis no assumptions have been made regarding either applicable 
regions or the initial properties of the two interacting waves. Equation 3.10 offers a 
means of expressing the modulated short waves as a combination of spectral 
coniponents based on monochromatic input. The short wave is viewed now as a 
carricr subjected to a simultancous amplitude and phase modulation, (Panter, 1965). 


Changing the limits of summation, equation 3.10 becomes 


n' = a'Jg(m,)sin p. (3.11) 


3, (m,.) {sin(yotny,) + (-1)"sin(po-nyy)] 
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The first term has the form of the solution for waves at the incident wave frequency 
propagating on a sloping bottom but with a different Bessel argument (see cquaton 
2.16 ). The remaining terms represent the side bunds at sum and difference frequencies 
(We nWy ) Where they are phase locked. This solution is dependent on the urgument 
my, and converges rapidly as n increases. 

In shallow water, the Bessel function argument my, reduces to 2 (ay /h) k, x 
for the case of progressive short waves superposed on long standing waves. The 
arguinent my 1s found to increase onshore (lig 3.1), where the computed values of my, 
arc evaluated until closer to the shore (x = I cm) before my, reaches © . Due to tlie 
onshore increase of m, , J, (m, ) decreases onshore, oscillating rapidly and reaches 
zero at the shoreline (lig 3.2). Note the diflerenee in behaviour of the Bessel functions 
in equation 3.11 and 2.15, where the Bessel argument allows the amplitude of the long 
standing wave to reach a maximum at the shorelinc. For breaking short waves on a 
long period standing wave, the depth limited breakers is utilized to describe a, (= ¥ h) 


in equation 3.5. Substituting a’ inta 3.11, yiclds 


; 3 2 
n' = yh[Jg(m)sin py. + 7 pH I, (m,) 008 WC] 


; Nn. ! 
+ yh{J (m){sin(p. +mp,) - (-1) sin), -m),)} (3.12) 


a 
- 5 Hm) Thy (MH) {cos ty, +My) - (-1)"cos () =n, )] 

Two distinct features of this solution are demonstrated in the amplitude squared 
spectra (lig 3.3). Tirst, symmictric side bands are generitcd about the scarier (iige3.3 
a). The second feature is the transfer of short wave energy to the upper and lower side 
bands as the depth decreases (Mig 3.3 b), which causes broadening of the wave spectra. 
It 1s also noted that the amplitude of the spectral components decrease shoreward. 
The mstabilitics associated with the side bands generation are attributed to their spaual 
structure and phase variation. It is shown that the side bands energy is modified by 
the decreasing depth. The modulation may contribute to the infragravity waves by 
Strengthening the low frequency energy component closer to the shoreline, which has 
previously been observed. Waddell (1973) presents evidence of nonlinear transfer of 
energy from high to low frequency due to collision, where successive bores begin to 
overrun cach other. [fowever, Huntley and Bowen (1975b) suggested that energy 


transfers are dué to the interaction of many irregular spaced breakers in the surf zone. 
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Figure 3.1 Wavenumber modulation factor m, vs non-dimensional distance X/Xp » 
Where X) 1s the surf Zone width. 
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Sawaragr and Iwata (1974) studied experimentally the wave transformation after 
breaking and observed energy transfer from the carrier into higher frequency waves by 
an unknown mechanism. The present analysis may be used to illustrate some of these 
findings. Equation 3.1] shows that the beach does not eliminate the incident wave 
frequency and that the modulated wave energy components and the incident wave 


frequency vary spatially. 


B. CHANGES IN THE WAVE KINEMATICS 


1. Background on Wave Kinematics 
At any instant of time, a wave front 1s defined by the phase function y = 
constant, which is the equation of a fanuly of parallel planes with normal vector k. As 
time increases, these planes move with the phase speed in the direction K. 


In general, the dispersion of tlie wave motion may be written 


w = w(k,h) (3.13) 


where the frequency varies with both the wavenumber and the local water depth. The 
local wavenumber and frequency are defined by the gradient of the slowly varying 


phase function in space and time: 


(3.14) 


xt 
T 

<] 

< 


wo = - 3 (outs) 


Since the curl of a gradient is zero, it follows immediately from equation 3.14 that the 


local Wavenumber vector in space 1s irrotational, L.c., 
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Figure 3.2. Behaviour of the Bessel functions J 
(solid curve) and J, (Mp) (dashed curve) across the 
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one =250 (3.16) 


(3.17) 


Equation 3.17 is known as the kinematical conservation of wave density, where the rate 
of change of wavenumber ts balanced by the convergence of the frequency. Sometimes 
3.17 ts also called the “conservation of crests equation,” since the crests are neither 
created nor destroyed and their total number must be conserved. 

Rewriting equation 3.17 in tensor notation and expanding the convergence 


term using the dispersion relationship in equation 3.13 gives 











a : kk : = ) 
ai, ou , eee h, Ro eee ee (3.18) 
gt ge OX 5 oh OX 5 


Equation 3.18 can be rewritten as 
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Figure 3.3. Amplitude square spectra of the modulated short wave 
at offshore distance: (ajx = 40.0 m(b)x = 0.5 m. 
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2. Kinematics of waves on current 
In studying the problem of short gravity waves superposed on waves of a 
much longer wavelength and period, the long wave may be modeled as a slowly varying 
current U(x,t), (e.g. Longuet-lliggins and Stewart, 1960, 1961; Garrett and Smith, 
1976). A coordinate system that moves with the current velocity is chosen, by which 


the observed frequency @ , passing a f1xed point, 1s given by 


a = We een (3.21) 


Where @, is the intrinsic frequency, (Ursell, 1960; Whitham, 1960). The subscripts s 
and I refer to short and long wave properties. The second term on the right hand side 
of 3.21 is,due to the relative notion “ORsioltumaves Withmenespectete the chosen 
coordinate and is known as the doppler shift effect. 

The short wave will be modulated by the presence of the current, and 


therefore the modulated short wave (denoted by primes) 1s assumed to satisfy equation 
Saleee. 





= +Vu' = 0 (3.22) 


where k’ , the observed (modulated) wavenumber, is unknown and assumed unsteady. 


The observed frequency is given by equation 3.21. Expanding equation 3.22 gives 





ye (3.23) 


Where the quantities k; , cg and @ are for the short waves, and U ts the current (long 
Waves). 


Assuming a steady intrinsic short wave train; 1e., 0k, / Ot = 0, equation 3.23 
may be simplified by the use of equation 3.19 to give 
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ee = 6 (3.24) 
x 


which describes the kinematic conservation of the modulated wave density. 


3. Applications 
The following analyses demonstrate the efficiency and simplicity of equation 
3.24 in examining the wavenumber changes of short waves riding on a current of 
longer period. Two different bottom topographies are considered: flat bottom and 


plane sloping bottom. Different current models are applied. 


a. Waves on a Current in Water of Constant Depth 
P 


The current will be modeled as a progressive or standing long period wave. 


1. Current as a progressive long wave 
The short wave is assumed to propagate over a flat bottom and to be 
parallel to the long wave (collinear) in the X-direction. The current U 1s described by 


the orbital velocity of a progressive long period wave, 


eosin k, (h+2) 


z ee ee oS ene ea! Pes. 
Ue “e_ “sinh k,h Q (3.25) 
Which corresponds to a free surface clevation, Ny = ap sin Wo Substituting into 


equation 3.24 gives 





Aen cosh k,(h+z) ; =e 
* +k f[ajw,k, — SS Sin b,)] = 
ae SAR Sinh kh A (3.26) 
Equation 3.26 is easily integrated over time to give 
ee cosh ky(h+z) ¥, = Constant anon 
s Sse sinh kh ae eeet 
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where the constant 1s evaluated at U = QO, so that k’= constant = k,. ,1.e. the 


S bf 
wavenumber will not change if there 1s no current. Then, the general solution ts 


obtamed such that 


S Cc, Ca28) 


which describes the wavenumber changes due to the non-linear interaction between the 
short wave and a progressive long wave using a kinematic approach. [Evaluating 


equation 3.25 at the free surface, z = 0, equation 3.28 may be written 
kK’ = ko ( 1+ ag kpcoth keh sin Wy ) 


[t is interesting that the above result, obtained from the first principles of kinematical 
conservation, 1s the same equation as (3.4 ) obtained by Longuet-Iliggins and Stewart 


(1960) using a lengthy second order perturbation analysis scheme. 


il. Current as a long standing wave 


For the flat bottom case, equation 3.24 reduces to 








The current U 1s described by the orbital velocity of a linear standing wave 


cosh k, (h+z2) 


mee ee t 
0 2 oe an kph Se (3.30) 


which varies only in the x-direction. Applying equations 3.29 for both x and y 


directions yields 


akY cosh k, (h+z) . 
se * Ky 82%, —ginh kph COS KyX Sin mgt = 0 (3.31) 
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Integrating Over time gives 


cosh kK, (h+2) 


eg ame k —_—_— — 
x S Ko Spot teen np Constant 
a : (3.32) 
Ke = k&k = €oOmstanc 
Be y 


Equation 3.32 indicates that the wavenumber component alongshore is always constant 
since no longshore component of the current exists. 
The constant in (3.32) can be evaluated at any node whiere Np = 0. At 


these locations, it is assumed that the wavenumber does not change, i.e. ky = k 


S 
« . e . « « ° X 
Substituting into equation 3.32, a general solution is obtained: 
Re = K . {l + a,k,coth k yh cos kx cos wt] , 
“ (3.33) 
fom). Kk ="BConstant —. 
M Sy 


Therefore the wavenumber component in the x direction is always modulated in time 

and space, while the longshore wavenumber component is always constant. As a 

result, the wavenumber vector changes in time and space; 1c. the waves are refracted. 
It is clear that if the waves are parallel to the current, 1c. collinear in 


x-direction, then k,, = 0 and the solution reduces to the first equation in 3.33 . 


cs 


b. Waves on Current in Water of Variable Depth 
A reasonable model to describe the current on a sloping bottom ts the 
orbital velocity of either a cross-shore long standing wave or edge waves. Collinear 
and oblique incident waves on a current are included. Let us assume a plane sloping 


bottom for simplicity of analysis. 


1. Current as a cross-shore long standing wave 

On a sloping bottom, the horizontal velocity of a standing wave 1s 
described by equation 2.17 . Recalling the equation for wave density conservation 
(3.24), 
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ak 
dU a 0 
Soy "a : + U - = 
aot =e Ox S ox 


(Sme=)) 
which requires specifying @ U/d x and dk, /@ x, . The gradient of the current velocity 
becomes 


an 
aU _ ae oh A =) 
Aeon x ie) eee 


components gives 


ey) 
Using the dispersion relationship in shallow water to differentiate the wavenumber 


dK. F , 
x oe 
a OT Ae 


and dk, / 0 y 


(3.36) 
Q, since the dispersion relationship is independent of y as the depth 
contours are assuined straight and parallel. Applying Snell’s law to describe the 
change in the angle of approach with x, d @ / @ x gives 

_ = 
ely 


oh tan o@vean sf 


Substituting into 3.36, gives 





(3.37) 
ok 
SX _ Sx Z 
Teens, SSS (3.38) 
Subsututing into 3.34 using 3.35 and 3.38 
ak! an k Scene 
a + k [- U - el By - U[ SX 
ot Sy x he ot 


(3.39) 


3k! 
7 ns 


, =k, = constant. Integrating cquation 


from which ¢ ky /d y = 0 and, therefore, k’ y 


So ONer UMC PIVCs 





Rex seca Ne 
ro ~V a ; = ~~ l= 
Ky = UU, tebe |[ tL ee aaa } Lae h Constant (3.40) 
since 
f U dt = (1/a? yd U/ dt + constant 
Therefore, equation 3.40 becomes 
k n 
: sx 3 Seea si man se 
NS = 5 aq li ts 5 ] a h Constan (3.41) 
XU) 9 


The constant in 3.41 can be evaluated at xX = © where current vanishes, ic. Ns = 0 
and U = 0, which gives k’, = constant = k, . Then, for a stcady oblique wave 
X 


train, the gencral solution is given by 





Z 
n 
_— a SCS td 
x ue 
[f the waves and the current are collinear (@ = 0), equation 3.42 reduccs to 
‘ 1  9U 
ee, ft 5 a ae! 
| x Kid 5 (3.43) 
i eons tan& 
4 ¥ 


This equation shows the dependence of the wavenumber changes on the surface 


elevation and on the acceleration of the current. 
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ii, Current as an Edge Wave 
In the 3-1) model the long wave is described by an edge wave. A zcro 
mode edge wave is used to simplify the analysis where the Laguerre polynomial equals 


unity. Phe horizontal velocity components are given as 


a,gk. ~k X 
U = ees e MY Simik - yu. ©) (3.44) 
W, yY Q 
Me x 
5G > x 
— dF yt d =r: e 
V 77 e Sea Wo ) (3.45) 


A steady short wave train riding on an edge wave will satisfy the wave density 
conservation equation (3.24 ). Substituting into equation 3.24, using cquations 3.38 


and 3.39, gives 


k 
ok! S 
——= U[k k ate = Becom = 0 (3.46) 
at ex YQ 
aky 0 
st Se eek (3.47) 


Integrating both 3.46 and 3.47 with respect to time gives 


‘: soem (3.48) 
Kk) + ko (eS = eons 
x Sy, Wy 2Xuw 7 


The constant in 3.48, can be evaluated when x =» ©O where both U and V go to zero 


and Noe = he . This leads to 
bre 


4$ 


2 
_ | a See 3.49 
MS Kk, tt (ls See (3.49) 
x R He 
y Y 
<= ee 
pe ec om (3.50) 


Since edge waves are progressive in the alongshore direction, equation 3.50 is similar to 
3.28 where the current was modeled as a progressive wave. The changes in the 
wavenumber and direction of short waves riding on zero mode edge waves are 
demonstrated in Pig 3.4. Vhe modulation in the wavenumber increases onshore 
Where the wavelength goes to zero. Increasing the angle of wave incidence 1s shown 


to increase the Wavenumber modulation. 


C. CHANGES IN WAVE ENERGY AND MOMENTUM FLUX 
1. Background 
A wave train composed of a group of individual waves propagates with the 
group velocity at which the total energy ts transmitted. Energy 1s usually expressed in 
terms of avcrage energy over a complete wavelength per unit surface arca. The total 
energy of any wave system is the sum of both the potential and Kinetic energy. 


The average potential energy L per unit surface area is given by 


iL Ne eee lB > 
Sy eae) ! 5 
BE, = 3 n J a ie (3.51) 
The average kinetic energy E, per unit surface arca is estimated by 
Ee, ele n=0 
2 | f (uve 4w")dz dxdt E32) 
ae 2Tk., . a 


lor linear progressive surface gravity waves, the average total energy [oy 15 
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Figure 3.4 Spatial changes in the wavenumber of shert period waves 
due to the presence of zero mode edge wave. 
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1 
En = 50g a (3.53) 


which 1s equally partitioned between the potential and the kinetic energy. 


a. First Order Sea-Swell Energy Modulation 
Short period waves riding on much longer wave can be exainined using the 
principle of superposition. First order modulations are examined using the linear 
description of both waves. The linear theory is recognized as a_ rcasonable 
approximation for its successful use in describing many observed phenomena even 
within the surf zone, (Bowen et al., 1968; Guza and Thornton, 1980). The total surface 


elevation may then be expressed as 
= -+- ; 
Nep ee (3.54) 


Where the subscripts s and 1 refer to sea-swell waves and long period waves 
(infragravity Waves) respectively. Inside the surf zone, short waves break, dissipating 
their energy onshore. The depth-controlled wave breaking model is adopted (equation 
2.11) to describe first order short wave dissipation with zero amplitude at the 
intersection of the mean water level (MJWL) with the beach. The depth in equation 
2.11 is modified to include the slowly varying depth changes due to the presence of the 


long wave, 


2 an oa Bian (3.55) 


where Np Is the surface elevation of low frequency wave motion on a sloping beach and 
is described by either a long standing wave or an edge wave. The objective of this 
approach is to examine a first order energy perturbation (second order in amphtude) 
due to the presence of long waves. In the 2-D model, the long wave ts modeled as a 


long standing wave. The potential energy for the combined short and long waves 1s 
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calculated by substituting the combined surface elevation given by equation 3.54 into 
equation 3.51 . The kinetic energy (equation 3.52) is similarly calculated using the 


total horizontal velocity, 


un = UL tu, (3.56) 
where 
v= y(h + ny) Yg/h cos wt (eo7} 


and u, is given by 2.17. It is assumed that the vertical velocities in shallow water are 
small compared with the horizontal velocities, i.c., w << u. The total energy L-p 
averaged over short and long wave periods is the sum of both the potential and kinctic 


energy and 1s given by 


1 2 1 eee oe 
G) x, 
(3.58) 
+ aes 36.00 | 
(Wo +W,) 

232 
toy a. ol (x) | ] 

4 ek WW, 


The total averaged energy in equation 3.58 1s represented by four terms at frequency 
bands @, ,a) ,(M,. + @)) and (@, - a) ). The first term in right hand side is identified 
as the short wave contribution if the short waves were considered alone. Similarly, the 
second term represents the contribution from the long waves alone. The third term 1s 
due to the nonhnear interaction between the short wave and the long wave and is 


represented by two side bands at the sum and difference between short and long wave 


frequencies. 
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b. Model Parameters and Preliminary Results 

Energy calculations for the combined short and long waves requires 
specifying the two parameters Y and a) in equation 3.58 . The breaker index y , which 
relates the breaking wave amplitude to the local water depth, is taken to be .20 when 
Psi sGerOOt Wein sqaure wave hem (Thorton and Guza, 1982)) The second 
parameter, the long standing wave amplitude, is defined as the vertical excursion of the 
low frequency swash oscillation at the shoreline. Larlier work by Munk (1949) and 
Tucker (1950) showed the infragravity wave height to be 10 percent of the incident 
wave height, with both quantities measured in 1S m water depth. Goda (1975) 
measured infragravity wave heights in 1 m water depth and found them to be 20-40 
percent of the offshore incident wave height. Guza and Thornton (1982) measured 
infragravity run-up height in Torrey Pines to be 70 percent of the wind wave heights 
measured in 10 m water depth. Holman and Bowen (1984) measured the significant 
runup height, defined as the significant vertical excursion of water level at the 
shoreline, and found it to be 60 percent of the incident significant wave height at the 
breakerline. Guza and Thornton (1985a) analyzed data acquired from three different 
experiments and found that the significant vertical swash excursion on the average was 


equal to the incident significant wind wave height measured in 10 m depth 


fe= Ho. (3.59) 


Therefore, the amplitude of the long wave ag is taken as .50 Hig . This result will be 
used here since equation 3.59 1s weighted over three different beaches and incorporates 
the Santa Barbara data, which are considered in the present analysis. 

Data acquired at the Leadbetter Beach experiment in Santa Barbara during 
[950 aremminizeceto specify the parameters of the energy model. The best {it to the 
beach profile assuming a plane sloping beach is found to be 1:25 with a surf zone width 
of SO m. The characteristic periods for both the short and long waves as determined 
from the spectral analysis of current meter records are 14 and 85 seconds respectively. 
These are typical values for the data acquired on I'ebruary 4 (narrow banded waves). 
The model results are shown in big 3.5. The potential energy for long waves 1s shown 


to reach a maximum at the shoreline with a magnitude comparable to the maximum 
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short wave potential energy at the breakerline (Fig 3.5 a). The kinetic energy for both 
short and long waves decays shoreward (Fig 3.5 b). The total averaged energy decays 
onshore, with short waves dissipating all their energy at the shoreline where the long 
waves have their maximum energy (Fig3.5 c). The two side band components show a 


steady growth (order of 15 percent) in total energy in the vicinity of the shoreline. 


2. Mass and Momentuin Fluxes 
a. Background 

The changes in the mass and the momentum fluxes are required to specify 
the driving forces for predicting the surf zone dynamics. Conservation equations 
developed by Phillips (1966) are followed, in which the terms representing the mean 
and fluctuating quantites have been separated. These equations are applicable to 
wave Inotion as well as general turbulent motion. No restrictions are placed on the 
wave slopes and/or amplitudes through the use of these equations. 


The conservation of total mass per unit areca can be expressed 


oM. 


30D i | eye 3.60 





Where 1,) refer to horizontal coordinates, D 1s the total averaged depth of water, which 


may include the wave sctup, and M, is the total mass flux defined as 
i saa * 2 
M = i 6) us Za, el ’ (3.61) 


Where u; is the total horizontal velocity that contains a mean quantity and a 
fluctuating quantuty due to the waves. The overbar denotes , at this stage, averaging in 
time over the short wave periods, but not so long us to exclude the long waves. 

The momentum flux equations are derived from the horizontal momentum 
cquauions integrated over depth, averaged over time and using the free surface and 


bottom boundury conditions: 
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aM. 





ao ee i. 2 Oe (3.62) 
one ss a OS aa! oe OX; a 


The first term on the right hand side is the change in the mean hydrostatic pressure 
due to the slope of the mean free water surface, and RR; is the mean shear stress term. 
The first term on the left hand side is the local change in horizontal momentum flux; 
the second term contains the mean momentum flux U; Mi and the excess momentum 


duc to the wave presence S;. , which was termed the “radiation stress” tensor by 


Longuct-Iliggins and suceeded 1962). The advantage of using the “radiation 
stress” technique to solve physical problems is that the second order effects are 
obtained using first order wave theory, and therefore some non-linear wave propertics, 
such as setup and surf beat, can be explained casily. The radiation stress for a 


progressive linear waves is given by 





- _1_} 
13 oD (3.63) 


where 6 is the Kronecker delta. The third term can be neglected since it is of higher 
order in deep water and has a negligible contribution in shallow water. Evaluating all 


terms in equation 3.63, the excess momentum tensor for a linear progressive wave is 


Be kK. 1 Cg . 
ij = =! ee ty Ee ag 16; Pomme int lene (3.64) 


where Kk, Ky are the wavenumber components and are used in the expression to 
mndicate the directions k, /k = cos @ and ky /k = sing. 

These expressions are also a good first approximation for slowly varying 
depth if a and k are interpreted as the local values corresponding to the local depth h 
(Mei, 1983). In shallow water, the group velocity Co equals the phase speed ¢ and the 


radiation stress components are approximated by 
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(3.65) 
2 rue oe 
oy = a sin a (3.66) 
S =e Same COS 0 (2.00) 
yx 


As expected, for a standing wave on a flat bottom, the radiation stress components are 


exactly twice the values of a progressive wave (Longuet-lIliggins and Stewart, 1960). 


b. Changes in The Radiation Stresses Due to the Long Waves 

The objective now is to rederive the excess momentum due to combined 
short and long waves to develop explicit formulae. The short waves are assumed to 
arrive obliquely at an angle @ , measured counter-clockwise between the wave ray and 
the positive X-axis. The long waves are modeled as a normally incident standing wave 
on a plane sloping bottom. ‘The procedure of calculating the momentum flux 
components due to the presence of both the short and long waves ts to substitute into 
the general formulation 3.63, which is applicable to all kinds of steady and unsteady 
flow. In the following development, the terms second order in amplitude (first order 
im energy and momentum) are retained, and all the higher order terms are neglected. 
Then, the resulting expressions are averaged over the long wave period to elunimate the 
time dependance of the obtained formulae. 


Phe radiation stress component Shs defined as the flux of onshore directed 


momentum, ts given from cyuation 3.63 as 


[ae ol 1 - 
ght J (puu+P)dz - 5 pgD (3.68) 


The first term under the integral is the product of the total horizontal wave velocities, 


which can be obtained by summing equations 2.6 and 2.17. The limits of integration 


2 


are taken from the bottom to the sull water level (SWL) at 1) = 0, which is consistent 
with the analysis being second order in wave slope. Integrating over depth and 


averaging in time over long wave periods yiclds 


va) 


2 2 1 


© ho 


0 i Z 
= —: h 
f puu dz xg 1 
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The pressure term in equation 3.68 is obtained by integrating the vertical momentum 
cquation (neglecting the vertical stresscs) from any depth z to the free surface (sce 


Phillips, 1966), which gives 


NY A n 5 n 5 
P | = of f gdz+ = f wdz + — f uwdz - ow” } (3.70) 
Zz eZ 


The first term is the hydrostatic pressure. Integrating from the bottom to the free 


surface and averaging over short and long wave periods gives 


NO 
No 





2 
1 2 ieey Die 
term l = pg a =i (nth)* + ws if. 1a, J, (x) } 3.71) 


The second term ts similarly evaluated, but it was found to vanish when averaged over 


the long wave period. The third and fourth terms give 


iL a es 2a 
term 3 = = qoha,w, Ig (x) -J, (K)] (3.72) 
1 w,tan 8 vgn Bey (eee 
+ Go abies g if 0 A iL ak 
(O75) 
22a Zoe 
7 pee ae aeFote re 2 owes 
term 4 = eae ae ( 4 5 Wo 5 


The time averaged radiation stress Syy 1s obtained by adding cquatrons 
Bao9. 3.71, 3.72 and 3.73: 


a, 2 Ze 
ogt3y7h? + (1 + *hyatagGo + Zaps) 


a (3.74) 


le 2 A | aay Se ; a oe 2 3 
+ qPewga, (tan 68 vYgh JQ (x) Ty (X)) oe oe h 


2 
Gs) 
Dp 2 XR 2 vs Z i Z 2 - 
4p sa hia (dg (x) +07 OO) = Zwodg (x) } 


The first term on the right hand side is identified as the short wave contribution if 
considered alone, which is in agreement with cquation 3.64 obtained by 
Longuct-Iliggins and Stewart (1962). The second and the third terms are the 
contributions by the long waves and the side bands. 

The transverse component of radiation stress, Sux, is defined as thie onshore 
flux of longshore momentum across a vertical plane any distance offshore. This 
component can be calculated from the direct evaluation of the momentum flux as the 


product of the total horizontal velocities, up and v7, 
Sn = RACE (3.75) 


A cross-shore long wave is considered for simplicity, i.c. no longshore velocity 


component, While an oblique short wave is considered. Phen, equation 3.75 becomes 


eee oe 
8. = 9 J coe Cosa tu, ) us Suquiel (al yz (3.76) 


The transverse component S,, 1s the driving force for the longshore current inside the 


yx 
surf zone. When considering the cross-shore long waves, the primary change in Syx 


arises from the amplitude modulation of the short waves, which is a function of the 


2, 


long waves. Integrating equation 3.76 over depth and averaging over the short wave 


period gives the time dependent form, 


iL 2 2 . 
S = segly (h +n) Isineae cos 


yx Cr 


Averaging over the long wave period to climinate the time dependency, the transverse 


component S,, inside the surf zone 1s 


VX 


_" 1 2) a 
oo = sogty he ee sy a,Jo(x) }sin a cOS a 


(3.78) 


The first term on the right hand side is in agreement with equation 3.67, while the 
second term can be identified as side bands. In equation 3.78 , it is noted that the 
cross-shore long wave by itself does not contribute to Syx , but their interaction with 


the short wave results in the side bands that increase the Syx term. 
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IV. DATA ANALYSIS 


The analysis of field data is performed to test the hypotheses that the nearshore 
dynamics are duc to both the short and long waves. Outside the surf zone, waves 
arriving In groups cause a change in MWL (Fig. 4.1 ) and drive a forced long wave 
(Longuet-Higgins and Stewart, 1962). At the breaker zone, the short waves break and 
their heights become depth dependent, obeying the saturation curve. The water depth 
is modificd by the presence of the long waves. Therefore, the short waves are expected 
to be modulated by the infragravity waves inside the surf zone. The experiments are 


briefly discussed, followed by analyzed results. 


A. DESCRIPTION OF THE FIELD SITES 

Two experiments were performed as a part of the Nearshore Sediment Transport 
BitdyeiNols), The objective of NSIS was to develop an improved engineering 
formula to predict sediment transport on beaches. Freld sites were selected for their 
simple topography of straight and parallel depth contours, which match the 
assumptions often used for nearshore dynamics. The first experiment was located at 
Torrey Pines Beach (T.P.), San Diego, California and conducted for one month in 
November, 1978. Torrey Pines is a gently sloping beach (tan Pp = 0.002) with nearly 
straight and parallel contours. The bottom is composed of moderately sorted, 
fine-grain sand (mean diameter 0.15-0.2 mm). The waves break as spilling or muxed 
spilling-plunging. A wide variety of wave and weather conditions were encountered 
during the experiment, from small to large (2 m) waves, from very narrow-band swell 
to wide-band sea, from calm to windy days (> 10 m/s). A total of 42 sensors were 
deployed over an area 520 m long parallel to the shore and up to 500 m offshore. 
Surface elevation and horizontal, orthogonal velocity components were measured along 
a shore-normal transect from offshore at the 10 m depth contour to the shore. 

Leadbetter Beach, Santa Barbara (S.B.), California was the second NSTS 
experimental site (Iebruary, 1980). Leadbetter Beach has relatively straight and 
parallel nearshore depth contours with well-sorted fine to medium size sand. The mean 
nearshore slope varicd between 0.017 and 0.05 during the experiment, depending on the 


wave climate. Well-developed cusps occurred during the beginning of the experiment. 
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Figure 4.1 Schematic diagram of short and long wave models outside 
and inside the surf zone. 
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No major offshore bar system was apparent. The shorcline bas an unusual cast-west 
orientation along a predominantly north-south coast. The open ocean waves are 
limited to a narrow window of approach (+ 9° centered on 249°) because of 
protection from Point Conception to the north and the Channel Islands to the south. 
The generally highly filtered ocean swell type waves from almost duc west must make a 
right angle turn to approach the beach normally. The result is often a narrow band (in 
frequency and direction) spectrum of waves approaching at large oblique angles to the 
bottom contours in the surf zone and driving a strong longshore current. Rod and 
level surveys of the beach profile were obtained daily. An array of current meters and 
pressure sensors were deployed extensively in a similar manner as at T.P. Severe 
storms occurred during the experiment, causing significant changes in the shoreline 
configuration and beach profile. Despite the destruction of most of the instruments 
during the storm, a very substantial and valuable data was recorded (Gable, 1980). 
Data were acquired at a sample rate of 64 Hz for several hours each day, then 
low-pass filtered and reduced to 2 Hz. The sensors, data acquisition system. and 


experimental sites are fully described by Gable (1979, 1980). 


B. DATA ANALYSIS 

Data acquired from T.P. and S.B. have been extensively analyzed by different 
investigators. Incident wave analyses include shoaling, setup, swash oscillation, 
runup, velocity moments, longshore current and sediment transport.  Infragravity 
waves Were studicd for runup, identifying both the cross-shore and longshore structures 
and type of the dominant motion, whether they are leaky and/or edge waves (Guza and 
Thornton, 1982,1985a,b; Oltman-Shay and Guza, 1986). As mentioned before, the 
dynamics of the infragravity waves were not previously studied. Furthermore, their 
effect on the incident waves, especially inside the surf zone, is of interest to understand 
the nearshore processes. I[t is hypothesized in this study that the interaction between 
the short and infragravity waves gencrates other waves (side bands) at the sum and 
difference frequencies. These waves occur at the neigbourhood of the short wave 
frequency, and they are of importance in driving longshore current and sediment 
transport at the shorcline. The existence of both the infragravity waves and the side 
bands are investigated utilizing records of wave staffs, current meters and pressure 


sensors. Coupling between the incident waves and the infragravity waves are tested via 
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spectral and cross-correlation techniques, and different theories are discussed. Since 
the present developed model is monochromatic, only days of narrow banded waves are 
considered. ‘Vhree days are selected from cach experiment (total of 6 days): 3, 4, and 6 
February 1980 at S.B. and 10, 20, and 21 November [978 at T.P. 


1. Spectral Analysis 

Velocity and pressure time serics of 68-minute Iength are uscd to compute 
spectra with 32 degrees of freedom and a resolution of 0.004 IIz. The power spectra 
are computed for all the available working sensors at cach selected day. The spectra 
are plotted on linear scale to demonstrate the relative tmportance of cach peak and 
normalized by the variance (area under spectra). In Fig 4.2 , a well defined swell peak 
is seen ata frequency around f, = .06 Iz in most of the instrument locations during 
10 Nov., T.P., and in Figures (4.3, 4.4) at freguency f, ~ .07 IIz in all the analyzed 
instruments tn 3 and 4 Ieb., S.B. A consistent valley near 0.05 Iz 1s present in almost 
all spectra, separating the incident wave band from the low frequency band. Within 
the tnfragravity band, a pronounced low frequency peak ts shown in most of the 
spectra at I = O17 [lz for 4 Feb,S.B. and .0195 Iizm 10 Nov.,1T.P. 

‘the mfragravity energy has been studied extensively since the first 
observations by Munk (1949) and Tucker (1950). Iluntley et al. (1981) identified a 
progressive low mode edge wave tm the longshore direction that satisfies the dispersion 
relationship of edge waves. Oltman-Shay and Guza (1986) extended the work of 
Huntley ct al. (1981) and were able to determine the energy contained in each edge 
wave mode. The low mode edge waves (n S 2) were found to doniunate the longshore 
current energy While leaky mode or high mode cdge waves dominated the cross-shore 
direction. They attributed the daily changes tn infragravity cnergy levels to the 
variation m meident wave conditions and the resonant forcing of edge waves described 
by Gallagher (1971). Guza and Thornton (1985a) observed surf beat in 3 field sites 
including $.B. and T.P. They found a phase shift of m /2 between the surface clevation 
and the horizontal velocity within low frequency motions (f S .O5 Iz) which mdicates 
the predominance of the standing wave structure of surf beat. This result agrees with 
Suhayda (1974). Guza and Vhornton (1985a) rejected the hypothesis that surf beat ts 
primarily generated by bore-bore capture within the surf zone (as hypothesized by 
Bradshaw, 1980 and others) because bore-bore capture models do not predict incoming 
surf beat energy outside the surf zone and do not describe the gradual tncrease in the 


infragravity band onshore. 
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Figure 4.2.) Normalized power spectra at different locations in [0 Nov., TP. 
showing an onshore relative increase in the infragravity energy level. 
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The transformation of the incident wave spectra across the surf zone ts shown 
in Figures 4.2 and 4.3 . A dramatic shoreward tncrease in the Iow frequency band 
energy is accompanied with rapid attenuation in the progressive tncident wave energy 
duc to breaking as thre depth decreases. The infragravity energy can exceed the wind 
wave energy by a factor of 2:5 im Wetnner stimiezoneie, metic a), oc) ne stonna 
conditions, Ifolman et al. (1978) observed a considerable increase tn the infragravity 
energy levels responding to the changes in incident wave conditions. It is clear from 
Figures (4.2, 4.3 and 4.4) that the spectra closer to the shoreline 1s dominated by the 
low frequency band. 

The theory of wave-wave interaction presented in Chapter 3, by which side 
band waves are generated as a result of the non-linear interaction, can be tested using 
the available data. Llgar and Guza (1985), on studying the nonlinear dynamics of the 
shoaling waves using bispectral analysis, showed that the low frequency motion at the 
beach face is significantly non-lincarly coupled to higher frequency modes seaward the 
surf zone around the peak of the power spectra. They did not resolve possible side 
band gencration in their analysis, choosing instead to obtain higher stability (low 
resolution). Using higher resolution here, peaks at frequencies (f, + ) are found in 
Fig. 4.5 which are identified as upper and lower side bands. This supports the 
hypothesis of the short wave modulation derived in equation 3.11 . Evidence of side 
bands is present for almost all the measurements on all days analyzed, but usually only 
the upper side band at frequency (fp + {9 ) is present (see Iigures 4.2, 4.3 and 4.4). 
Vhe side bands are not always statistically significant, however, since at least the upper 
side band is present for all narrow band days. Ilenece, tt ts concluded that the side 
bands are real. It is not clear why the lower side band ts absent. The suppression of 
the side band may be strong lower frequency breaking or it may be caused by the 
seaward radiation of energy by the outgoing component of the standing long waves. 
LiIgar (1985) showed that the low frequency reflected waves will have phase 
relationships relative to their incoming parts that tend to decrease bispectral level. 
Guza ct al. (1984) show an inercase mn coherence between the incident wave envelope 
and the long standing waves when removing the reflected wave. 

The time history of energy transferred between swell peak and the upper side 
band is demonstrated in Fig. 4.6 using almost a four-hour time series. The long record 
is divided into 64 segments cach of 17.1 mmutes with a 50% overlap to calculate the 


spectra. I:nergy appears to be transferred from the incident wave peak to higher 
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frequency components. This leads to a faster growth of the upper side band than for 
the lower side band. The incident wave peak drops to a minimum, strengthening the 
upper side band. A similar mechanism could happen if lower side bands exist, which 
may transfer energy to low frequency components. 

Using the same data set, Guza and Thornton (1985b) computed the velocity 
variance in 20 Nov., T. IP. at each instrument location across the surf zone Fig. 4.7 . 
Their result show no substantial differences between the observed and the predicted 
values When summing the potential and kinetie energy for short waves riding on long 
standing waves (Fig. 3.5). The decrease in short wave energy is almost balanced by the 
increase in low frequency energy resulting in an almost constant cross-shore velocity 
variance (Guza and Thornton, 1985b). Wright et al. (1982) obtained a similar variance 
distribution on Australian beach. Although the spectral shape of 3 Feb is similar to 4 
Feb, waves in the 4 Feb are more energetic (Hein = 80 cm). Therefore, strong 
non-linearities are expected. 

2. Cross-Correlation Analysis 

The objective of this analvsis is to cross-correlate the incident wave envelope 
with the long waves to test the hypotheses that long waves are forced outside the surf 
zone and short waves modulated inside the surf zone. The cross-correlation function, 
Cp , is a measure of the degree of linear relationship (association) between two data 
sects in the time domain. Tucker (1950) correlated the incoming wave envelope and the 
corresponding long wave observed by a bottom mounted pressure sensor at a distance 
800 m offshore. He obtamed a negative maximuin correlation corresponding to a 5 
minute time lag. The time lag associated with the maximum correlation was 
approximately the time required for the incident wave group to reach the shore with 
, and for the associated long wave to travel back with phase velocity // gh. 


g 
The negative correlation was described by Longuct-Higgins and Stewart (1962) as a 


velocity c 


low frequency forced wave response due to the incident wave groupness that has a 
phase shift of m with the envelope of the high frequency waves. Ilence, under high 
waves, a long wave trough 1s expected and vice-versa. 

Iluntly and Kim (1954) obtained a positive correlation between long waves 
and the short wave envelope at zero time lag from current meter measurements taken 
outside the breaker linc. Their results agree with the foreed wave model of 
Longuct-Higgins and Stewart (1962). Tor velocity measurements, a positive correlation 


is expected with the x-direction positive offshore. 
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Figure 4.6 Spectral time serics of current meter CODx with 50% overlapping 
calculated from four hour record, (4 Feb, 1980, Santa Barbara). 
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Guza, ct al. (1984) obtained a negative correlation outside the surf zone and a 
positive correlation inside the surf zone when comparing the surface elevation of the 
short wave cnvclope with long waves. No explanation was given for change in signs of 
the correlution between outside and inside the surf zone. 

Surface clevatron time serics are inferred from current meter data from 
experiments at T.P. and S.B. using the complex Fourier spectra of the horizontal 
velocity component, U(f). Current meter records were first linearly detrended to 
exclude the effects of the rising and falling of the tides, then high-pass filtered using a 
cut-oll frequency of 0.05 IIz to exclude surf beat, followed by a low-pass filter with a 
high frequency cutoff (0.311z) to exclude capillary waves. The filtering was 
accomplished by [lourier transforming the signals and zeroing the Fourier amplitude 
coefficients in the filtered frequencies. The complex surface elevation spectrum, X(f), 


is calculated applying the lincar wave theory transfer function, I¢(f) 


(eee H,(£)U(f) 
(4.1) 


Then, the complex surface elevation spectrum 1s inverse transformed to obtain the 
surface elevation time series. The entire 68-minute-record was transformed at one trme 
to minimize the end effects, which result in spectral leakage, and to obtain maximum 
resolution for very sharp roll-off at the filter cut-offs. Surface elevations were also 
computed from pressure signals in a similar manner by transforming the pressure 
records using linear theory. 

Using the surface clevatron trme series obtarned in this manner, the correlation 
between the incident wave envelope and the corresponding low frequency motion 1s 
calculated. The low frequency motion is obtained by low-pass filtering the surface 
elevation time series. The high frequency wave signal is the residual after subtracting 
off the low frequency energy. The short wave envelope is obtained by demeaning and 
squaring the high frequency signal, followed by the same low-pass filter used before. 
Segments of the time serics for the wave envelope and the assocrated low frequency 
motion suggest positive correlation for the wave sensor tmside the surf zone (I*ig. 4.8) 
and negative correlation for the current meter derived signal outside the surf zone (Fig. 
4.9). 


The cross correlation function is calculated using the following formula 


12 


300 400 500 600 


depth (cm) 


200 





© © © © 
= & 8 & 
( ,208/,Wd)< o> 


f elocity variance in 20 Nov., 1978. 
Low is (0<1<0.05 Tz) and high is (0.05 <f<.5 Elz) 


Figure 4.7 Spatial distribution of band passed v 
| (after Guza and Thornton, 1985b). 
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where tT = k tis the time lag, t is the sampling time interval, n is 
number of data points, 6) GO, arc the standard deviations for the low frequency signal 
and the envelope, and x, y are the means for each of the time scrics. The time lag is 
positive when the long waves are in advance of the wave cnvclope, and it is negative 
when long waves lag the wave envelope. 

A test of significance (null hypothesis) based on the “fFisher-Z 
Transformation” is performed with the statistic z = Bes) 2 he (lt i l- Cp) 
(Miller and Freund, 1977). At level of significance @ = 0.05, the calculated z always 
exceeds 2,,/2 (= 1.96), since the used record is long ($192 data points). Thus, the null 
hypothesis is rejected and the correlation is significantly different than zero. 

The correlation function for the wave staff W38 (T.P.) ts shown in Fig. 4.10. 
The maximum correlation is +0.49, corresponding to a time lag of -4.0 sec. The 
function shows a strong periodicity at about 40 sec. The correlation coefficient 1s 
calculated with 95%. confidence to be significantly different from zero. Tig. 4.11 shows 
a negative correlation for CO7x (S.B.), indicating the characteristics of a forced wave 
response outside the surf Zone. 

The maximum correlation coefficients and the corresponding time lags at cach 
instrument location in TO and 21 Novemiber, |.P. are given in lable J and 11 
Similarly, tables 11] and 1V are given for 4 and 6 Feb.,S.B. Lower values of Cp inply 
that the two random variables may be independent, that there 1s no linear relationship 
between the two signals, suggesting a free long wave (Iluntly and Kim, 1984) or the 
possibility of some non-linear relationship. The low values of the negative Cp can be 
attributed to the fact that the m phase shift is only a deep water case, while the deepest 
Instrument used is 3m. The other reason could be the significance of the phase-locked 
reflected low frequency waVes decreasing the scolicrencesc cn ior low valucs ol 
reflecuion coefficient (LElear, 1985). Llvarwand = Guz (1955) Sinai tecmen oremeipe 
infragravity modes do not appear to be bound with a fixed phase relationship to the 
high frequency wave groups, since their biphase evolves towards lower values as waves 
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Figure 4.10 The cross-correlation function for W38, 21 Nov., 1978, 
Torrey Pines. A positive maximum correlation (.49) corresponds to -4.0 time lag. 
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Iigure 4.11 The cross-correlation function for CO7x, 4 Feb., 1980, 
Santa Burbara. A negative maximum correlation (.29) corresponds to + 3.0 time lag. 
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It is of interest to observe a consistent pattern of Cr , being negative outside 
the surf zone and positive closer to the shoreline icine 2 sarcdieetii3 jae | lie 
Statistically significant values of Cp in the nearshore zone at zero time lag indicated 
significant coupling between the high and low frequency motions. A zero correlation 


in Fig. 4.12 is found to coincide with the approximate location of the breaker linc. The 
zero correlation couklt be due to the two signals being non-linearity related, 


randomness (noise) introduced into the wave signal by the breaking wave, a possibility 
that the breakpoint could be at a nodal point for the long standing wave, independence 
of the two signals, or a transition location between forced and free long waves. Fig. 
4.13 shows a consistent negative correlation outside the surf zone and a positive 
correlation shoreward of the midsurf zone. This could mean that the wave breaking 
does not destroy completely the wave groupness until some location mside the surf 
zone. This location nmught be a transition between the forced and free long waves. 
Elgar (1985) showed that surf beat biphases decreases as waves shoal. On some days, 
Meewcorrclation in Pig. 4.13 shows a rapid decrease when Waves get closer to the 
shoreline. This might be attributed to the bore-bore capture mechanism. Mase and 
[wagaki (19S4) illustrated that a wave (bore) could be captured by a successive bore 
before the bore reaches a maximum runup height. Therefore, the number of runup 
waves and frequency of runup will be reduced compared to that of incident waves. The 
rauo of number of waves ts given using the surf sinularity parameter and shows a 
@eercase in number of waves reaching a gentle slope beach. 

The positive correlations obtained tnside the surf zone are demonstrated by 
the simple model developed in chapter 3 by linear superpose the two waves. The 


surface elevation time Serics mside the surf zone is simulated using 


n =) heeeme ia 


i )cos Va am (4.3) 


L 
The wave envelope and low frequency motion are obtained as described above. A 
positive correlation is always shown at zero time lag (fig. 4.14), which agrees with the 
observations. 

In the following chapter, the theoretical results obtained in Chapter 3 are 
applied to derive wave setup, longshore currents (steady and unsteady) and sediment 


transport due to the combined short and long period waves. 
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Figure 4.12 Spatial distribution of cross-correlation coefficients 
showing a dominant negative correlation outside the surf zone | 
and a positive correlation inside the suil cone. B. L. 1s the treakerline. 
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Figure 4.13 Spatial distribution of cross-correlation coefficients 
showing a dominant negative correlation outside the surf zone 
and a positive Correlation after midsurf zone. 
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Iigure 4.14 
betwecn the wave envelop 


V. APPLICATIONS 


A. WAVE SETUP AND SETDOWN 
1. Background 

In the nearshore zone, it is observed that waves during shoaling and breaking 
produce a variation in the mean water level (MWL). This variation in MWL is 
considered a primary cause of nearshore circulation, such as rip currents. The change 
in MWL, 1, is characterized by two regions (a) a gradual depression of the sea water 
level starting from offshore, and reaching a maximum at the breakerline, and (b) an 
increase in the MWL shoreward of the breaker line. The depression is termed 
“setdown” and the upward slope ts called “setup.” 

To study the wave setdown and setup, the X-component of the horizontal 
momentum flux equation 1s considered where the waves are assumcd to be steady state 


and propagating normal to the shoreline over straight and parallel depth contours, 


a 


9S. -, on 
a See ee (say) 





Equation 5.1! indicates that the change of the excess momentum {lux due to wave 
motion in the cross-shore direction is mainly balanced by the mean water slope dn/dx, 
assuming no gradients in the Y-direction and neglecting the frictional forces. 

[t is common to separate the nearshore zone tnto two dynamical regions of 


“outside” and “inside” the surfzone. 


@ieOutsidesthe Surf Zone : 

[f waves are assumed to propagate onshore without any energy loss, then the 
onshore component of the radiation stress 1s expected to increase steadily as the depth 
decreases. Therefore, the mean water level 1s lowered (sect down) by the presence of the 
unbroken waves. Longuet-Iliggins and Stewart (1962) integrated equation 5.1 by 


assuming that y < < h to obtain 


See eee (5.2) 
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This equation 1s second order in local wave amplitude. The setdown increases steadily 
until the breaking zone. Without referring to the radiation stress, Longuct-I liggins 
(1967) also derived the above solution using Bernoulli’s integral equation. This simple 


technique will be illustrated and applied below. 


(2) Inside the Sumi Zone 
As waves break, the wave amplitude inside the surf zone is assumed to be 
controlled by the total local water depth including the wave setup yn , simulating a 


spilling breaker type, 1.c. 
Ay eee (58) 


The cross-shore component of the radiation stress, S is a function of the energy 


XX? 
density which decreases inside the surf zone and reduces to 3/2 E in shallow waters 


SUC Ciabe 


ae 2 — 2 
qogy (h +n) (5.4) 


—N 


Ax 


Both the gradients of S,, and the forcing in equation 5.1 have a negative slope. 
Therefore, the change in the MWL inside the surf zone is expected to be positive, 


representing a wave setup. Integrating equation 5.1 gives 


ee SS) ¢ Ny (5.5) 


where N = (1+ 2/3 ye yy! . The subscript b indicates the value at the breakerline 
where Nh =- y Il, /8 is determined from equation 5.2 . Equation 5.5 shows an 
increase in the mean water level with decreasing depth, balancing the decreasing 


radiation stress duc to wave dissipation inside the surf zone. It should be noted that 
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Bernoulli's integral outside the surf zone is not applicable inside the surf zone since the 
mowuon 1s no longer irrotational. 

Bowen et al. (1968) demonstrated the above theoretical results by careful 
measurements in the laboratory. Experimental results agreed well with the theory. 
Ilowever, just outside the breaking point, their experimental valucs were less than the 
theoretical predicted since the theory does not consider energy dissipation. Also, very 
close to the shoreline, Bowen ct al. (1968) found a residual wave height (possibly due 
to a standing Wave) to persist and corrected the above cquation by letting the sctup 
slope be tangential to the beach slope, ic. dn/dx —» - dh/dx as (h + yn) 0. Van 
Dorn (1976) found that the setup gradient at the shoreline is actually steeper than 
predicted by equation 5.5, which agrees with Bowen ct al. (1968). Battjes (1974) and 
Bates and Janssen (1978) used a random wave model, which predicts smaller setup 
values. 

In addition to the wave setup laboratory investigations, a few field 
measurements have been reported. Dorrestein (1961) measured the change in mean 
water level across the surf zone during times when the significant offshore wave height 
ranged fron: 0.8 to 1.6 m. Ile measured a maximum setup of 0.15 m. Guza and 
Thornton (1981) measured the setup at the shoreline when the significant offshore 
wave height varied from 0.6-1.6 m. The maximum sctup at the shorelrme was found to 
be about 0.17 Ig , where Il, is the deep water significant wave height. Ilolman and 
Sallenger (1985) measured the setup under conditions of incident significant wave 
height varying from 0.4- 4.0 m. They found the setup to vary linearly with the surf 
srmilarity parameter Cy = B (IIo /Lo y Hie 


2. Changes in Setup and Setdown Due to the Presence of Long Waves 
Changes in the wave sctup and setdown due to the presence of long waves can 
be demonstrated by considering the same wave ficld used before, 1.c. oblique short 
waves riding on a normally incident long period standing wave. The wave sctdown 1s 
evaluated applying Bernoulh’s integral equation (Longuct-Iliggins, 1967). Furst, the 


vertical momentum cquation 1s integrated over depth and averaged over ume. 


‘3 GnGia)s ee (5.6) 
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where p is the pressure and the overbar denotes mean quantities. Next, the time 


averaged Bernoulli's integral equation is evaluated at z = 0, assuming the flow is 
irrotational, 
P | a: 2 0G ve ae + a9 + Constant = 0 (5.7) 
Z—0 : z=0 


- & —( Gaeee—a + Constant 
n 7G i. (5.8) 


The difference in mean sea level 6 W is obtained by. applying equation 5.8 at two 


different locations (xX, yy ,0) and (X> yz ,0), which avoids evaluating the constant. 
sq = alto - at | (5.9) 
g 


For the combined short and long standing waves, equation 5.9 1s used to evaluate the 
setdown as follows. Each velocity component is expressed as the sum of the short and 
long standing wave velocities as given before in chapter 2. Let location 2 be in deep 
water (1) —» © ) where the setdown is zero. Substituting into equation 5.9 and 


averaging over the long wave period gives 


2 
eS alee, 2 


- pat Se eee (5.10) 
2 sinh 2kn  2g‘“g ~ Mg 


5 | 
i 


where subscripts s and | refer to short and long waves. The first term in equation 5.10 
is identified as the the setdown solution (equation 5.2) when considering the short 
wave only. At the’ breakerline, the depth controlled breaking model is utilized and 


setdown reaches a maximum to give 
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where the subscript b refers to conditions at the breakerline. 

Wave setup can be evaluated by spatially integrating the X-momentum flux 
equation, where the radiation stress component S,, for the combined wave fields is 
given by equation 3.74. The constant of integration is evaluated by applying the 
obtained setup equation at the breakerline where the change of the mean sea level 1s 
known from equation 5.11. The resulting equation is quadratic in 4 and has two 
roots. One root is disregarded since it unreasonably predicts setdown at the shoreline. 


The second root gives 


ee oih | Yh -O ) (5.12) 


where 
aie — u Jinnae a 2 
= - [ny sy) hong + 51 Sy (hy ) 
, (5.13) 
is 
eC Oe 3, (XY) 
,° 2 0. 
ig 
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Using Longuet-lliggins’ setup model for short waves as a reference, the setup is 
relatively increased closer to the shoreline due to the presence of the long standing 
wave (Fig 5.1 ). Ilowever, the setdown is also increased at the breakerline compared 
with the reference model. The relative increase in setup appears to be over-predicted 
and suggests the importance of including damping terms, such as frictional dissipation 


and/or percolation, in the swash zone. 


91 





0°08 O°0Z 0°09 O°0S OOF ODE O'0Z O'OT OD 
(WI) dALAS SAUM 


7 5.1 Wave setup due to both incident and infragravity waves 
compared with honeuer venta and Stewart model (1962). 
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B. LONGSHORE CURRENT 
|. Background 
When waves obliquely approach the surf zone, a mean current parallel to the 
shoreline is generated, which is important in transporting sediment alongshore. In the 
simplest formulations, the longshore current is assumed steady and two dimensional, 
being mdependent of the longshore topographic changes. The linear Y-momentum flux 


equation 1s simplified such that the gradient of the radiation stress component @ Syx {Ox 


1s balanced by the bottom shear Stress, 


6 (5.14) 





Outside the surf zone, the wave energy is assumed non-dissipative and the excess 
momentum flux in the Y-direction is a conservative quantity. Therefore, no driving 
force 1s expected outside the surf zone for generating a longshore current (Bowen, 
1969). Inside the surf zone, cnergy is dissipated due to wave breaking, causing a 
change in the momentum flux. The onshore changes in momentum tflux drive a 
longshore current primarily confined to the surf zone. 

Nuimerous models of longshore current on a long straight beach have been 
done since the introduction of the radiation stress concept. For obliquely incident 
waves on plane beaches, analytical solutions were derived for monochromatic waves by 
Bowen (1969) and Longuet-!Iliggins (1970 a,b) and for random waves by Thornton and 
Guza (1986). Numerical treatments (e.g. Thornton, 1970; Jonsson et al., 1974; 
Madsen et al., 1978); Wu and Liu, 1984, and others) have employed arbitrary beach 
profile and/or including setup and/or non-linear inertial terms. In monochromatic 
wave models, lateral mixing is introduced to smooth out the discontinuity in the 
longshore distribution due to the intensive current shear at the breakerline (Bowen, 
1969, Vhornton, 1970; Longuet-Iliggins, 1970a).. For random waves models, there 1s 
no need to include the lateral mixing (Collins, 1970; Battjes, 1972) since the waves 
transform and break over a considerable horizontal distance resulting 1n a smooth 
change tn the rate of energy dissipation and Syx 

Periodic variation and unsteadiness in the longshore current measurements 
have been reported by Wood and Meadows (1975) and Meadows (1976), but neither 


U3 


analytical nor numerical solutions have been attempted. Mcadows (1976) suggested 
three distinct longshore velocity components, a steady component and two fluctuating 
components, with one at a short wave frequency and the other at low frequency. 
Guza and ‘Thornton (1978) and Holman and Bowen (1984) observed low frequency 
components in the longshore velocity spectra. Since the longshore current is a 
time-averaged velocity, Guza and Thornton (1978) stated that an appropriate temporal 
averaging time for mean longshore current is unknown. The above studics suggest that 
the low frequency wave motion can be important inside the surf zone. 

Infragravity waves are included in the following wave description tnside the 
surf zone, allowing a more complete description of the wave field. The changes tn Syx 
due to obliqueness of the short wave and the total bed shear stress are derived. Then, 
the steady and unsteady longshore current formulations are considered on plane 


sloping beaches. 


2. Wave Refraction and Radiation Stress 
The depth contours are assumed to be straight and parallel, and obltquely 


incident waves are therefore refracted according to Snell's law, 


= eee = eee (Sl) 


where c 1s the phase speed. The S,, component for the combined short and long 


yX 
waves is given by equation 3.77. Assuming a small angle of wave tncidence (cos @ ~ 
1), Snell’s constant (sin « /c), which contains all the wave angle information, can be 


introduced into equation 3.77 to yield 


(5.16) 


Another approach for deriving equation 5.16 can be shown by using equation 3.67, 


where 
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BE 
S = —+ sin a cos a 


yx (5.17) 


By assuming a small angle of wave incidence, 5.17 is reduced to E Co (sin a /c) where E 
is given by 1/2 p g y2 (h+ ny ) . This will give the same solution as equation 5.16. 
This result can be attributed to the fact that standing waves by themselves do not 
contribute directly to the longshore momentum flux since they are oriented in the 


cross-shore direction, and no S,.. forcing is expected for a normal standing wave. The 


va 
contribution of the infragravity Waves manifests itself by the generation of momentum 


at the side bands. Averaging equation 5.16 over the long wave period gives 


S Bs ete an fo OC + 5a JQ X)) (5.18) 


yx 2 Cc = 


That is to say, the S,, forcing is increased by the formation of the side bands (second 


term in bracket). The longshore current forcing term 0 Syx / 6 x for both steady and 
unsteady cases 1s compared in I‘ig. 5.2 with the Longuet-I liggins (1970a) modcl where 
only short waves are considered. In general, the combined forcing by short and long 
waves increases in the vicinity of the shoreline due to the relative growth of the side 


bands. 


3. Bottom Shear Stress 
Inside the surf zone the bottom stress can no longer be neglected, especially 
when considering the change of momentum flux in the Y-direction. Assuming the 
quadratic shear stress formula, the averaged alongshore bottom shear stress, exerted by 


the oscillating wave motion Is given by 





<= 2 aac (5.19) 


<< 


where cr is the friction coefficient and uy is the vector sum of the current and the 


inviscid orbital velocities for both short and long waves measured just above the 
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Figure 5.2 Gradient of the longshore current forcing across the 


Surf zone. 
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bottom boundary layer. The absolute value is necessary to insure that the stress 
reverses with velocity reversals. Longuet-[liggins (1970a) assumed a small angle of 
wave incidence and a weak current compared to the water particle motion of the 


waves, and he was able to simplify 5.19 to 


a a “ 3.20 
7 31. [Py V (5.20) 


where V is the mean motion in the longshore direction which is assumed uniform over 
depth. The weak current assumption has been shown to be a _ reasonable 
approximation for many field applications (Liu and Dalrymple, 1978). The total 
velocity vector uy can be resolved into its components in the cross-shore and 


alongshore directions: 


ue Sie a) as (5.21) 
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Where u, and uy are the X-component of the orbital velocities and 1,j are the unit 
vectors in the x and y directions respectively. U and V are mean velocities. Assuming 
uniform currents alongshore, the mean cross-shore velocity U is zero according to the 
conservation of mass, since the beach forms a boundary in the X-direction. 
Considering an oblique short wave riding on a cross-shore long standing wave, 1.¢. Vy 


= Urequatioms.21 can be rewritten 





cette cosa + u.)}i + {(u. sing) + wea (S722) 
oh S ye S 

The absolute value of the total velocity vector is 

| U.. | = Vu Cosa + uy)” fee S Pico ) /) (5.23) 


| 


Assuming a small angle of wave incidence (@ < 10° ), cos @ can be approximated to be 


unity, Which simplifies 5.23 to 


2 
ml = Cu. te Uy) + ve + 2u_v Sin a (5.24) 


= (uv 4S) eel (5.23) 


which can be rearranged into 





ele al (40 (5.26) 


Note that uy + uw = uy 1 according to equation 3.56. Then, using the weak current 


assumption, V/uy < <1, to simplify the analysis gives 


(5.27) 


Since the two velocities u, and wu, are almost collinear (both acting almost in the 


X-direction), equation 5.27 gives the basic relationship 


— i ste rea] (5.28) 


Equation 5.28 can be evaluated using 3.57 and 2.17 and averaged over the short wave 


period to give 


| u = /g7nt Li + a5 J a eal +a, |Jy x) sin w,t|} 


7 (5.29) 
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The alongshore bottom shear stress averaged over the short wave period is then given 
By the Usero! >. 20", 


ee 
tT = pe, Vg/hi=*(|h+ a,Jp(x)| +a, 15, (x) sin wet|}V (5.30) 


V 1 


Since 5.29 is still time dependent, the expression is averaged over the long wave period 
to fit the steady formulation 


uy, = eV eh (yh + aio +5, OH. (5.31) 


The corresponding longshore bottom stress is given as 


Z eT cag . . oe 
ee oc. O/ Ril yh+ ay (y¥Jg lx) | +]Iy Cx) |iV (ha) 


4. Longshore Current Models 


Obliquely incident short waves riding on a normally incident long standing 


wave are considered to develop steady and unsteady longshore models. 


a. Unsteady Model 


Neglecting the nonlinear mean advective terms, the linearized Y-momentum 


flux equation has the form 


n 0S 
JM bs ne ee iy (333) 
oe 0 


where My = V p his the mean longshore mass transport averaged over the short wave 


period. No variation in the longshore direction 1s allowed in equation 5.33, and the 


oe) 


lateral mixing 1s ignored for simplicity. Assuming an unsteady longshore current 
appears to be more general, especially with the presence of irregular oscillatory waves 
(Wood and Meadows, 1975). Thus, equation 5.33 can be written in terms of the 


longshore current V: 
valine (5.34) 


The onshore divergence of the radiation stress S,, is considered the significant forcing 


yx 
for generating longshore currents and 1s given by differentiating 3.77 , 


mye ey). = Weipa) + C™ (Nels ad (5.35) 
ame a 


k * * % ; * : 
where N is aconstant and B ’C andD are functions of space only. N_ is given as 


Va? Sauce aa 
ae (= eh 





N* = 
Cc 
(5.36) 
2 > tang nae 
and 
BhYGT (x) — 2H’? wd, (X) 
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The longshore current forcing, Syx /O x decreases onshore, then increases closer to 
the shoreline, having a singular point at x = O (sce Fig 5.2 ); the computed values 
stop atx = 2cm. The forcing term is shown to oscillate at the long wave frequency 
(fig 5.2). Dividing equation 5.34 by p h and calling the resulting forcing term q” 
(x,t), equation 5.34 can be simplified to 


SV g pe(x,t)v = g* (x,t) (5.37) 


(Oj) 2, 


Q? 


where p. Is given by 
3 ee . i 
pr(x,t) = ¢, Vg/n*t2t{h ta,gy00 | tadt, Od sin w tl? (5.38) 


It can be shown that the time dependent term (second term between brackets) is small 
compared with the short wave contribution in the swash zone, since uy / lu, | > 0 at 


the shorchinc. Equation 5.38 then reduces to 


Oe (sx, er Vg/n7 {24 |n + agg (x I) (5.39) 


Equation 5.37 is completed by specifying the initial condition V(x,0) = 0. An 
analytical solution to the unsteady longshore current is obtained using an integrating 
fectemeeMespite the fact that both the forcingetcrm and the bed shear stress have a 
singular point at the shoreline, the resulting longshore current 1s well behaved and has 


a finite valuc at the shoreline. The gencral solution can be expressed in four terms as 


wie, t) = Vsteady + Vunsteady l + Vinsteadv 2 (5.40) 


3 Veransient 


The first term is a steady term and 1s given by 


2 
meee, Sit) Oye ean es 2 
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where A” =y(Cht ay Jo (x)) . The steady term is similar to the 
steady solution given by Longuct-Iliggins (1970a) differing only by the second term in 
the A’ expression. The second term in 5.41 1s an additional steady component as a 
result of the side band formation. The two unsteady terms are functions of the 


frequency of the long wave and Its first harmonic and are given by 


ee” 
7 7 fae 5.9.9, (sin Oy 
Us ceady ol 2[4 cfg Are = 12 woh] e 
: (5.42) 
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The fourth term in equation 5.40 is the constant of integration 
and 1s evaluated by applying the initial condition at time zero. This term 1s identified 
as a transient term Which grows very rapidly to a steady state within a few long wave 


periods. 
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In the vicinity of the shoreline, the longshore current 1s found to be finite and steady 


and 1s given by 


tang Saeee 
Vo 
0) £ (5.45) 
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The spatial distribution of the longshore current is shown in Fig 5.3. The dominance 
of the side bands and the oscillating components starts at the most shoreward 1/5 of 
the surf zone width. At the breakerline, intensive shear occurs as expected in 
monochromatic models. Lateral mixing, if added, will diffuse the shear and smooth the 
longshore distribution as described before. Comparison with Longuet-lliggins’ (1970a) 
model shows an onshore increase of the longshore current and a relative decrease 
closer to the breakerline. The decrease is due to the the increase in the bottom shear 
Stress. 

The time variability of the longshore current is evident in Fig 5.4 , where 
the current is shown to oscillate at the infragravity wave frequency. In the vicinity of 


the shoreline, the first harmonic (V ) dominates the fundamental frequency 


unsteady2 
component. The transient term Is Just an initial response to the wave system that 


becomes neghigible within a few infragravity wave periods. 


b. Steady Model 

l’or steady flow, equation 5.14 is used in which the cross-shore gradient of 
S,., 1s balanced by the bottom stress. As in the analysis of the unsteady model, wave 
setup is not included. Steady solutions with and without considering lateral mixing are 
developed. Averaging over long wave periods is required to eliminate the time 
dependency in the equations describing the unsteady models. Vhe bed shear stress 1s 
given by equation 5.32 , and the transverse radiation stress component is obtained by 
averaging equation 5.35 over the long wave period. Substituting in equation 5.14, the 


steady longshore current 1s given by 


_ 2m sin Oy tanB 2 
steady 3 (e Cc eA** 
24? (5.46) 
Y 
7 ,sin a Q 2 
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Where A” = yh + a (y Jo (x )+ J, (X )). The first term of 5.46 1s similar to the 


. *% . 
steady term in the unsteady case except for A. At the shoreline, the growing side 
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bands will have a steady and finite contribution to the longshore current given by 
equation 5.45, t.c. for both steady and unsteady models the magnitude of the current is 
lic sammerat the shoreline. 

Ignoring the exchange of momentum duc to horizontal turbulent eddy 
transport increases the current shear at the breakerline and results in a sharp 
discontinuity m the velocity distribution. The presence of any horizontal mixing 
couples the adjacent water elements together and diffuses the momentum in the 
cross-shore direction. Adopting the concept of eddy viscosity, Longuct-Iliggins 


(1970b) expressed the mixing term to be proportional to the mean current V such that 


(5.47) 





where Vv = Np x ¥ gD is the horizontal eddy viscosity, D 1s the total water depth and 
N is non-dimensional coefficient, (0 < N < 0.16). Thus, the steady Y-momentum flux 


equation becomes 





mee, ty p24) (5.48) 
ax y ox ox 


Equation 5.48 can be rewritten as a second order differential equation tn V such that 


G Syy / @ eee a(x) 6? V/ o x2 + My (x) 6 V/d x + m3 (x) V (5.49) 
where 
m, = Nx hv gh 
my = SJ2NhV gh 
m3 = - pcr] Uy | 
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Tigure 5.4 Time variability of the longshore current 
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and |uy | is defined by equation 5.31 . Equation 5.49 is solved numerically using a 
finite difference scheme. The coefficients my, , My and m, are evaluated at mid-points. 
A tri-diagonal matrix is generated over a domain equal to twice the surf zone width. 
Two boundary conditions are specified at the shoreline and away from the breakerline. 
At the shoreline, the eddy viscosity term goes to zero, i.e. the mixing length goes to 
zero at a boundary, and equation 5.48 reduces to 5.14. Therefore, the steady solution 
without considering eddy viscosity, equation 5.45, is used to specify the boundary 
condition at the shoreline. This approach is necessary because equation 5.49 is 
singular at x = 0, and care must be exerted during the numerical formulation to avoid 
mathematical inconsistency. The longshore current is assumed to vanish at a distance 
of one surf zone width seaward the breakerline. 

Steady longshore current models induced by short and long waves are 
shown in Fig. 5.5 compared with Longuet-Higgins (1970a) model. A comparison with 
Longuet-Higgins (1970a) model shows a relative decrease in the current seaward the 
inner surf zone. This can be attributed to the increase in the total horizontal velocity 
associated with the bottom shear stress. It can be seen that including the lateral 
mixing smooths the longshore current profile at the breakerline and shifts the 
maximum velocity onshore. The longshore current distribution is similar to the results 
obtained from the EOF analysis by Guza et al. (1986) showing a non-vanishing current 


at the shoreline (Fig. 1.4). 


C. LONGSHORE SEDIMENT TRANSPORT 

Incident waves arriving obliquely to the breaking zone, releasing their energy and 
momentum, generate a longshore current which inturn drives the longshore sediment 
transport along the beach. Knowledge of sediment transport rate is essential to the 
understanding of various coastal engineering problems, such as topographic changes in 
beach profiles and the associated marine protective works. For estimating the 
longshore sediment transport, there are three general approaches: the wave power 
model, the traction approach and Bagnold’s energetics model. In the wave power 
model, the rate of total longshore transport Qp¢ is simply assumed to be proportional to 
the longshore component of wave power Po at the breakerline to some power n (Watts, 
LS 3), 
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Figure 5.5. Steady longshore Current models , with and without adding 
lateral mixing, compared with Longuet-lliggins model (1970a). 
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Oi el (5.50) 
where kK and n are empirical coefficients. Py is given by 


. 525) 
P, = (EC I) cosa, sin a), ( ) 


which can be expressed in terms of the radiation stress (Longuet-Higgins, 1972) 
P a C (5.52) 


Inman and Bagnold (1963) found n in equation 5.50 to be unity and suggested the use 


of the nnmersed weight sand transport rate 
(5.53) 


where Pp, , Py are the density of sand and water and a’ is the sand porosity (roughly 
0.60). Komar and Inman (1970) measured longshore sand transport rate, combined 
their data with other data to seck a general formula. Most of the data falls onto a 


single trend satisfying the simple relation 


a ae? Tee (ena) 


where Py» is in watts/m when specifying wave height in terms of TH... and Qe is in m? 


/day. Ilowever, Equation 5.54 1s empirical with no real consideration to sand transport 


mechanisms. [1 addition, the coefficient .77 has a high degree of uncertainty. 
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The traction approach includes work done by Biker (1971), Swart (1976), 
Madsen and Grant (1976) and others. For instance, Madsen and Grant (1976) adapted 
the Einstein-Brown sediment transport equations to time-varying sand movements 
under combined waves and unidirectional currents. The bottom stress and resulting 
sand transport vary with time. ‘heir model is relativly complicated and requires 
parameters for which no data are available. 

Bagnold’s energetic model assumes that the total sediment transport 1s the sum of 
the bed load and the suspended load transport. Bed load is maintained by 
grain-to-grain contact within 10 cm of the bottom while the suspended mode, 
containing generally finer grains, is lifted from the bed by the turbulent action of the 
wave-induced vertical motion. The suspended sediment grains are supported via 
turbulent diffusion. They move with nearly the local water velocity. Ilowever, they 
are also falling vertically relative to the local supporting fluid with a fall velocity, W. 
Thus, part of the available energy is extracted to maintain the particles in the water 
column. Bagnold (1963) developed the total load sediment equation for steady, 


two-dimensional stream flow, 


a i, + ig 
| _ E (5.55) 
= $$ 1 0) 
j tan ‘p tan 3 - or 3 | 


where 

1 = total immersed weight sediment transport (bed and suspended) 

fy = bed load efficiency 

fo = suspended load efficiency 

@ = internal angle of repose of the sediment 

tan f} = slope of the stream bed 

u = mean velocity of the stream 

W = fall velocity of the sediment 

(2 = rate of energy production of the stream 
Inside the surf zone (oscillatory flow), Bagnold (1963) relates the rate of immersed 
weight transport to the work done by waves and current on the bed surface. The 
mechanics of sand particle movement are described as a back-and-forth motion by the 


oscillatory wave motion close to the bed with essentially no net transport. Wave 
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energy is expended in supporting and suspending the moving sand above the bed. 
Once the sediments are suspended, the presence of any unidirectional current 
superimposed on the to-and-fro motion can then produce a net drift of sediment. 
Bagnold (1963) derived the following relationship 


ee (5.56) 


where ig is the time averaged immersed weight transport rate per unit width in the 
direction @ determined by the unidirectional current ug , £2 here is the available wave 
power due to energy losses by the fluid drag at the bed surface, u,,, is the magnitude of 
the orbital velocity of wave motion, and K’ is a dimensionless coefficient. ‘The success 
of the wave power and Bagnold’s energetics approaches is due to their relative 
simplicity, requiring only one or two parameters specified from the data. 

Equation 5.56 has been used as a basis for the development of a number of 
longshore transport models. Inman and Bagnold (1963) specified the various 
parametcrs of waves at the breakerline to calculate the longshore imimersed weight 


transport rate, 


=< 


K' (EC), COS ms) (5.57) 


rei 
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Equation 5.57 can be viewed as a general relationship because it does not specify the 
forcing for the longshore current. It could be due to a tidal current, a nearshore cell 
circulation, wind or to oblique wave approach. Komar and Inman (1970) required K’ 
to be equal to 0.28 to agree with their data. 

Bowen (1980) applied Bagnold’s model to the problem of on-offshore sediment 
transport on beaches with normally incident waves and onshore steady current (no 
longshore currents). Bowen used Stokes second order wave theory to predict the 
equilibrium beach profile as a function of incident wave characteristics. Using a similar 
approach, Bailard and Inman (1981) derived a general bed load equation in the 


lonyshore and cross-shore directions by considering time-varying flow over a sloping 
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bottom. Their solution reduces to equation 5.56 for steady, two-dimensional flow and 
assuming a weak current and a smiall ungle of incidence. Bailard (1981) extended the 
work of Bailard and Inman (1981) to include the suspended sediment load to give a 
more general expression. The total longshore sediment transport of Bailard (1981) 


solution is given by 


2 


2 
; ~ ee V 
_—— ~=—C} ean ae 
(5.58) 
€ 
S 3 
ieee: V 
ar wer W | | 


where the first term on the R.H.S. represents the bed load and the second term 1s the 
suspended load. The total average sediment transport can be obtained by integrating 


equation 5.58 across the surf zone from zero to the breakerline (Xx, ). 


ee emiinncs (5.59) 


Calibration of the model can be done using the wave power model. 

There is considerable uncertainty concerning the relative importance of bed load 
versus suspended load within the surf zone. According to field measurements, 
Komar(1978) suggested that the suspended sediment transport rate accounts for less 
than 20% of the total transport. llowever, Downing (1984) estimated the suspended 
load on a dissipative beach to account for as much as 45% of the littoral drift. 
Sternberg et al. (1984) measured the suspended sediment rate at Leadbetter Beach, 
Santa Barbara and found that the measured suspended sediment longshore transport 
rate is equal to the total longshore transport rate as predicted by equation 5.57. This 
result agrees with the data acquired by Kana and Ward (1981) during a storm. 
Sternberg et al. (1984) found maximum suspended sediment occured at the midsurf 
zone region, and pointed out that the suspended sediment appears to be high in the 
swash-backwash region. 

Kraus et al. (1981) measured the longshore sand transport distribution using 


multicolor fluorescent sand tracer; a distinct bimodal distribution was obtained. 
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Maxima in transport occured in the swash and/or within the breaker zone, where 
suspended load transport is expected to be relatively high duc to the turbulence. White 
and Inman (1986) measured the longshore sediment rate in Santa Barbara and also 
found a bimodal distribution across the surf zone with maxima near the break point 
and in the swash zone (fig 5.6). They estimated the longshore suspended load to be 
roughly 10-30% of the total longshore transport rate. 

The above results show a large volume of sediment in the swash zone which is 
not predicted by any available model. Waddell (1973) found a low-amplitude 
mugration of bed forms which migrated consistently downslope suggesting that uprush 
sediment transport was primarily through suspension sediment. Bradshaw (1980) 
pomted out that the swash-backwash collision process (low frequency waves) is 
significant and results in tenporary suspension of a large volume of sediment in the 
roll wave. Wright et al. (1982) demonstrated the dominance of the suspended sediment 
in the inner and swash zone using a turbidimeter situated 10 cm above the bed. The 
time serics shows pronounced turbidity maxima at infragravity intervals. In general, 
the predictions of littoral drift have a high degree of uncertainty, which results partly 


froin our basic inability to make accurate measurements of sand transport on beaches. 


!. Longshore Sediment Transport due to short and long waves 

Vhe incident waves are assuined to approach obliquely the surf zone at a small 
angle riding on a cross-shore long standing wave. [Tollowing Batlard (1981), equation 
5.58 is used to calculate the total longshore sediment transport distribution. Based on 
laboratory and field data, Batlard (1981) estimated the efficiency coefficients to be Ep 
= 0.13 and && = 0.025 and W =0.04 m/sec . The total horizontal velocity squared 
averaged over long wave period is used in evaluating the bed load term in equation 
oS 


C2) iin mais. aed bade) 


Similarly, the averaged cubic total horizontal velocity is evaluated to calculate the 


suspended load term, 
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The steady longshore current, with lateral mixing added, is 
utilized to derive the sediment alongshore. The longshore sediment transport 
distribution due to the combined waves is shown in Fig 5.7. Note that including the 
lateral mixing smooths the distributions of both the longshore current and sediment 
transport at the breakerline. The resulting distribution qualitatively simulates the 
bimodal distribution given by White and Inman (1986), since it describes the maximum 
sediments near both the shoreline and the breakerline. 
Unhke previous models, the wave energy has been shown to exist at the shore 
(in the swash zone), and hence observable dynamical processes can be described 
analytically. None of the available models for predicting the longshore sediment 
transport explain the increase in longshore sediment transport in the swash zone. All 
the previous models predict no transport at the shoreline which has been shown to be 
unrealistic. It is shown that the longshore sediment transport model presented here 
estimates correctly tle bimodal distribution and shows the dominance of the suspended 
load, which agrees with some of the above results. The main improvements of the 
presented model are the inclusion of the infragravity waves and, therefore, describing 


the swash dynamics and the merease of sediment transport at the shoreline. 


D. MODEL VALIDATION 

The wave-wave interaction approach within the inner surf zone suffers from our 
inability to model the breaking process accurately. Nonetheless it can lead to greater 
insight into the inner surf zone dynamics. Several assumptions are necessary for 
simphfying the analysis due to the complexity of real surf zone processes, but it 1s 
believed that they have some physical basis. It 1s assumed that the linear theory 
describes the waves adequately for the order of the analysis, and inside the surf zone 
spilling breakers are assumed, where wave height is strongly controlled by the local 
depth. Non-breaking long period waves on a plane sloping bottom are assumed to be 
either leaky mode or edge waves with a frequency 0) << @, . Basic assumptions of 
Straight and parallel beach contours have been used for simplicity where no longshore 


depth variation is allowed. Plane sloping bottom is assumed for simplicity. The 
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Figure 5.6 ‘Measured longshore sediment transport r 
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sensitivity of the model to the change tn beach slope 1s clear since infragravity waves 
modify the local water depth and, therefor modulate the total energy. Limitation to 
narrow banded incident and infragravity waves is required for the model input. It must 
be noted that using a monochromatic description of both incident and infragravity 
waves may raise questions about quantitative predictions since the processes in nature 
are stochastic. IIowever, the developed model gives a number of qualitative results 
which may be useful in understanding the inner surf zone dynamics. Application of the 
model to different nearshore processes leads to a simple physical picture of the nature 
of the wave-wave interactions in the inner surf zone. The validity of this model is 
confirmed by comparison with field data, where the resulting distributions at least 


qualitatively simulate actual distributions. 
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Figure 5.7 Longshore sediment transport distribution model due to both 
short and long waves showing a bimodal structure. 
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VI. SUMMARY AND CONCLUSIONS 


The modulation of incident short period wind waves superposed on infragravity 
waves are investigated. Modeling the short waves as depth limited breakers, the short 
waves are shown to be inodulated in amplitude, wavenumber and direction due to the 
slowly varying depth caused by the long standing waves. Furthermore, modulation 
results in additional energy at the side band frequencies near the incident wave peak. 
The energy calculations show a gradual growth of the side bands onshore with 
non-vanishing energy at the shoreline (order of 15%). An analytical expression is 
developed to deseribe the surface elevation of the modulated short waves as an infinite 
sum of spectral components at the incident and side bands frequencies. It ts shown 
that the modulated short wave transfers energy to the side bands as the depth 
decreases. Therefore, it is suggested that the spectral components of the lower side 
bands may strengthen the infragravity wave band as the depth decreases resulting in 
the infragravity waves dominating the power spectra. 

The changes in the wavenumber and direction of the short wave duc to the 
presence of a long wave are calculated by using a moving frame of reference that 
reduces the non-breaking long wave to a slowly varying current. On a plane sloping 
beach, the wavenumber modulation is shown to increase with increase angle of wave 
approach. The excess momentum flux components and bed shear stress duc to the 
combined waves are rederived analytically and show periodicity at the frequency of the 
long waves. 

Spectral analysis shows the dominance of the infragravity waves closer to the 
shoreline. Using narrow band data with high resolution (df = .004 I1z), side band 
peaks can be identified at frequencies (f, + f ). On almost all days and for all the 
instruments, the upper side band is well defined in the power spectra, but the lower side 
band 1s suppressed. The reason for the lack of the lower side bands on most days 1s 
not known. [fowever, breaking of the low frequency waves and/or radiation of energy 
by outgoing component of the standmg wave may cause suppression of the lower side 
bands. 

Most days show a positive correlation between the incident wave envelope and 


low frequency motron inside the surf zone, while a negative correlation is observed in 
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deeper water in agreement with Longuet-Higgins and Stewart's (1962) explanation of 
forced waves. The positive correlation is explained bv the short wave modulation due 
to the presence of long waves. However, the values of the measured cross correlation 
coefficients are found to be low (although statistically significant), partly due to the 
reflected wave field, that decreases the coherence (Elgar, 1985). On some days, the 
observed zero cross-correlation was found at approximately the breakerline and other 
days in the inner surf zone. This may suggest a nodal point for a standing wave and/or 
a transiuon location between forced and free long waves. 

It is hypothesized that some of the often observed discrepancies with present 
short wave dynamical models occurring in the inner surf zone are due to the interaction 
between the incident waves and the infragravity waves. The derived expressions for 
energy and momentum fluxes are uulized to study wave setup, longshore current and 
sediment transport in the inner surf zone. 

Wave setup calculations show an increase in the mean water level compared 
with setup due to the short wave only. This is because the breaking short waves are 
depth controlled and the depth 1s slowly varying due to the presence of the long Waves. 
Consequently, energy at side bands is generated and the setup is increased. 

A steady longshore current driven by oblique short waves riding on cross-shore 
long period waves, is found to be finite at the shoreline due to the presence of finite 
side band energy. The model, in general, predicts smaller currents than 
Longuet-Higgins’(1970a) model because of the relative increase in bed shear stress due 
to the presence of long waves. In the derived analytical expression, contributions by 
the short waves and side bands can be identified. It is shown that the cross-shore long 
wave does not contribute to the forcing term, 6 S.., / 0 x, since it is averaged out. 

It is more general to consider an unsteady longshore current. The analytical 
solution for an unsteady longshore current shows three major terms: a transient term 
that vanishes after a few long wave periods, a steady term similar to the one obtained 
in the steady model, and an unsteady term which oscillates at the long wave frequency, 
but driven by the momentum changes at the side band frequencies, and acts closer to 
the shore line. This unsteadiness is attributed to the side band formation. 

An improved longshore current model is necessary to improve the longshore 
sediment transport formulation. Based on the sediment transport model by Bailard 
(1981) and using the developed steady longshore current model, a longshore sediment 


transport distribution is calculated. The model predicts sediment transport maxima at 
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the shoreline and the breakerline. The results Qualitatively, the results compare 
favorably with various ficld data by Zenkovich (1960); Kraus et al., (1981), and White 
and Inman (1986). 

It is recommended to include the wave setup in the longshore current models. 
Numerical solution for the unsteady longshore current including the large angle of 
wave incidence and assuming a strong current is of interest. Numerical solutions can 
be compared with the obtained analytical results. Actual depth profiles are 
recommended to use in computing the longshore sediment transport rate to compare 


quantitatively with the data. 
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